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ri ! Abstract 



The Phenomenology Working Group at NOW'98 discussed the following topics: 

• Possible interpretations of neutrino oscillation data in the framework of 
neutrino mixing. 

• Phenomenological models that accommodate the neutrino oscillation data. 

• Ideas for future experiments. 

• The role of neutrinos in cosmology. 

• Quantum-mechanical problems of neutrino oscillations. 

• Problems of the statistical interpretation of the data. 
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1 Neutrino Mixing and Oscillations — Samoil 
Bilenky 

1.1 Introduction 

The neutrino mixing hypothesis |l|-0 (see also 0-0) is based on the assumption that 
neutrino masses are different from zero and the neutrino mass term does not conserve 
lepton numbers. In this case, the left-handed flavor neutrino fields u^l (a = e, /i, r) are 
superpositions of the left-handed components UkL {k = 1, . . . ,n) oi the fields of neutrinos 
with definite masses rrik: 

n 

T^aL = ^UakVkL, (1-1) 

fc=l 

where [/ is a unitary nxn mixing matrix. The left-handed flavor neutrino flelds v^l 
(a = e, /i, r) are determined by charged- current (CC) and neutral-current (NC) weak 
interactions with Lagrangians 

/:?^ = --^ Y, V^lpaLW + h-c, (1.2) 

/:r = - o \ E T^lpJ^aLZP + h.C.. (1.3) 

a=e,ii,T 

From the precise measurement of the invisible width of the decay of the Z-boson we 
know that the number of light flavor neutrinos is equal to three (see |]10|); corresponding to 



z/g, z/^ and Ur. On the other hand, the number n of massive neutrinos has no experimental 
constraint, besides being bigger or equal than three. 

There are two possibilities for the fundamental nature of massive neutrinos: 

Dirac. If the total lepton number 

L = Le + Lf, + Lr (1.4) 

is conserved because of the invariance of the Lagrangian under the global gauge 
transformation 

z/q ^ e**^ z/q , a^e^'^a (a = e,/i,r), (1.5) 

then massive neutrinos are Dirac particles. In this case: 

• The flelds u^ have four independent complex components. 

• It is natural to expect that the number n of massive neutrinos is equal to the 
number of flavor neutrinos, i.e. three, although nothing forbids in principle 
the existence of sterile Dirac neutrinos. 

• Neutrinoless double-/? decay {{PP)ou) is forbidden. 

• Dirac masses and mixing can be generated with the Higgs mechanism of the 
Standard Model. 



Majorana. If the Lagrangian is not invariant under the global gauge transformation 
( |1.5|) , the total lepton number L is not conserved and massive neutrinos are Majo- 
rana particles, i.e. truly neutral fermions which do not have any charge (electric, 
leptonic, etc.) that distinguishes particle from antiparticle. In this case: 

• The massive Majorana fields Vk satisfy the Majorana condition 

Vk = vl, (1.6) 

where z/^ = CVk^ and C is the charge-conjugation matrix. 

• Neutrinoless double-/? decay is allowed. 

• If right-handed neutrino fields VaR (singlets of SU(2)l) exist, the number n of 
massive Majorana neutrino is bigger than three. In this case to the mixing 
relations ( |1 . 1|) one must add the relations between the right-handed fields VaR 
and the massive fields Vk'- 



z/f 



aR 

k=l 



^Uak^kL- (1-7) 



The quanta of the right-handed fields are sterile neutrinos that do not partic- 
ipate to weak interactions. 

In the Majorana case there are two plausible options: 



(I) The see-saw option |^ . If the total lepton number is violated by the right-handed 
Majorana mass term at an energy scale much larger than the electroweak scale, 
the Majorana mass spectrum is composed by three light masses m^ [k = 1,2,3) 
and three very heavy masses Mk {k = 1,2,3) that characterize the scale of lepton 
number violation. In the simplest see-saw scenario (see, for example, |^,|T2|-|T3| and 
references therein) the light neutrino masses are given by 

^k--^^< (^ = 1,2,3). (1.8) 

where m^ is the mass of the charged lepton or up-quark in the k^^ generation. The 
see-saw mechanism provides a plausible explanation for the smallness of neutrino 
masses with respect to the masses of all other fundamental fermions. 

(II) The sterile neutrino option. If more than three Majorana mass terms are small, 
then there are light sterile neutrinos. In this case active neutrinos z/g, Vfj, and Ur 
can oscillate into sterile states. Notice that sterile neutrinos can be obtained in the 
framework of the see-saw mechanism with some additional assumptions ( "singular 
see-saw" [|I6|, "universal see-saw" W%)- 

The main open problems concerning neutrinos are: 

1. What are the values of 



(a) the neutrino masses mk7 

(b) the elements Uak of the mixing matrix? 

2. Which is the nature of massive neutrinos (Dirac or Majorana)? 

3. Which is the number of massive neutrinos? 

4. Can active neutrinos oscillate into sterile states? 

5. Is CP violated in the lepton sector? 

1.2 Neutrino Oscillations 

A neutrino of flavor a and momentum p produced in a weak interaction process is de- 
scribed by the state 

w.) = J2u:kWk). (1.9) 

k 

Here {uk) is a state describing a massive neutrino with momentum p and energy 

,2 



El, 



p^ + m^ ~ p + — ^ 



;i.io) 



where p = |p| . The expression ( |1.9| ) is based on the assumption that the state of a flavor 
neutrino is a coherent superposition of states of neutrinos with different masses. 

The general expression for the probability of z/^ — > z/^ transitions in vacuum can be 
written as 



l^a^l^0 



Sap + E ^"'^ ^/ 



f3k 



fc=2 



exp 



2E 



1 



;i.ii) 



where AmL = ml — m'j, L is the distance between the neutrino source and detector and 
i? ~ p is the neutrino energy. It is clear that neutrino oscillations can be observed only 
if there is at least one Am^j 



E 



> 1. 



;i.l2) 



There are three experimental indications in favor of neutrino oscillations coming from 
the results of 



1. Atmospheric neutrino experiments (Super-Kamiokande |T^, Kamiokande [19], 1MB 
Soudan p^) with the squared mass difference 



Ami„ ~ 3 X 10"^ eV^ . 



;i.l3) 



2. Solar neutrino experiments (Homestake [^, Kamiokande |2^, GALLEX p4 
SAGE [|^, Super-Kamiokande [^]) with 



AmL~10-'eV^ 



or 



Am ~ 

sun 



10-^^ eV' 



(MSW [27]) 



vacuum osc. 



(1.14) 
(1.15) 



3. The LSND experiment p8l, with 



Am 



LSND 



leV^ 



;i.i6) 



Furthermore, in order to extract from the experimental data information on the values 
of the neutrino masses and mixing angles it is necessary to take into account also the 
negative results of numerous reactor and accelerator short-baseline experiments (the latest 
and most restrictive ones are: Bugey ||29| for the v^ -^ ^e channel, CDHS [^ and 
CCFR [0 for the v^, -^ u^ channel, BNL E776 0, CCFR ^ and KARMEN g for 
the Vfj_ — i> Vf. channel, CHORUS [|35| and NOMAD |3^ for the v^ -^ v^ channel) and of 



the recent CHOOZ reactor long-baseline experiment [p7 



From Eqs. (|1.13|) - (|1.16|) the experimental results indicate the existence of three differ- 
ent scales of Am?, i. e. at least four massive neutrinos |3^^2| (see also Section ^ The four 
types of neutrino mass spectra that can accommodate the solar, atmospheric and LSND 
scales of Am? are shown in Fig. [1.1| . In all these spectra there are two groups of close 
masses separated by a gap of the order of 1 eV which provides the mass-squared difference 
AmLgNE, = Am\i = ml — ml that is relevant for the oscillations observed in the LSND ex- 
periment. Two years ago we have shown [3U| that, if also the negative results of numerous 
short-baseline neutrino oscillation experiments are taken into account, among the possible 



schemes shown in Fig. ^3] only the schemes A and B are compatible with the results of 
all neutrino oscillation experiments. In scheme A Am,l^^ = Arrigi = m,2 — mi\ is relevant 



for the explanation of the atmospheric neutrino anomaly and Am, 



2 
sun 



Amln = m? — m\ 



MS 



is relevant for the suppression of solar i/g's, whereas in scheme B Am 



atm 



Am43 ^^^ 



Ami^^ = Am,2i- These two schemes have important consequences for long-baseline ex- 
periments 1^, for the possibility to observe CP violation in the lepton sector |^ and 
for Big-Bang nucleosynthesis [H3,E3[. The phenomenology of neutrino oscillations in the 



schemes A and B is identical, but in the scheme A the effective neutrino masses in Tri- 
tium /3-decay experiments and in neutrinoless double-/? decay experiments can be of the 
order of m,^ ~ m^, whereas in the scheme B they are strongly suppressed ||39|. Hence, the 
observation of the effect of neutrino masses of the order of 0.1 — 1 eV in Tritium /3-decay 
experiments and in neutrinoless double-/? decay experiments could allow to distinguish 
between the two schemes, favoring the scheme A. 



1.3 Mixing of three neutrinos 

If the results of the LSND experiment will not be confirmed by future experiments, the 
most plausible scheme is the one with mixing of three massive neutrinos and a mass 
hierarchy: 

mi <C m2 <C ms . (1-17) 

In this case Amg^^ = Am|i = m^ — m,^ is relevant for the suppression of the flux of solar 
i^e's and Am^^j^ = Am|i = mg — mf is relevant for the atmospheric neutrino anomaly. 
The probability of Pa — * J^p transitions is given by 



^fp 



Yl ^»k Upk exp 



fc=i 



.Am|iL 
' 2E 



(1.18) 



-- 7714 



1712 

■^ mi 
(I) 



m4 

1713 
-- 777,2 



mi 
(11) 



7774 
7773 



.- 7772 
-- 777i 

(A) 



Figure 1.1 



7774 

77i3 



7772 
-^ TTii 



(B) 



Figure 1.1: The four types of neutrino mass spectra that can accommodate the solar, 
atmospheric and LSND scales of Am^. 



where E is the neutrino energy and L is the distance between the neutrino source and 
detector. Let us consider vacuum oscillations in atmospheric and long-baseline (LBL) 
neutrino oscillation experiments. Taking into account that in these experiments 



A7n2i-^ 
2E 



<1 



;i.i9) 



and using the unitarity of the mixing matrix, we obtain 

Sal3 + Uf33 ?7*3 



pLBL 



exp 



/\ml^L 



2E 



1 



;i.2o) 



Thus, under the condition (|1.19|) , the atmospheric and LBL transition probabilities in 
vacuum are determined only by the largest mass squared difference Am,'^^ and by the 
elements of the mixing matrix that connect flavor neutrinos with the heaviest neutrino 

From the expression ( |1.20| ) , for the probability of z/q, — > z/^ transitions with (3 ^ a and 
for the survival probability of z/q, we find 



P 



LBL 



*up 



jLBL 



1 ^ [^ AmliL 

2^.;. 1-cos^^ 



1 / /\ml.L 

l-2B„;Jl-C0.^|- 



for 



[3 ^ a 



with the oscillation amplitudes A^.q, and B^-a given by 

^/3;a = 4 I UjS^ I I t/a3 1 , 



;i.2i) 
;i.22) 



1.23) 
1.24) 



Hence, in the case of a hierarchy of the masses of three neutrinos, neutrino oscillations 
in atmospheric and LBL experiments are characterized by only one oscillation length. 




Figure |L| 
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Furthermore, from Eqs.( |1.2l| ) and ( |1.22D one can see that the dependence of the transition 
probabihties on the quantity Am1iL/2E has the same form as in the standard two- 
neutrino case. Let us stress, however, that the expressions ( [1.21| ) and (|1.22| ) describe 
transitions between all three flavor neutrinos. Notice also that the transition probabilities 
in atmospheric and LBL experiments do not depend on the possible CP-violating phase 
in the mixing matrix and we have 



P, 



LBL 



jLBL 



^up 



(1.25) 



As it is seen from Eqs.( |1.2l| )-(]TT^ 
in all channels (up ^ z/„, z/„ ^ u. 



in the scheme under consideration the oscillations 



'^M' ^f^ 



Z/p 



are described by three parameters: Am 



2 
31; 



It/p3p, |t//,3p (because of unitarity of the mixing matrix lUrsl"^ = 1 — It/^P — 1^,,'^'' 



With the help of Eqs.( |1.2^ ) and ( p..24| ), one can obtain bounds on the mixing param- 
eter |f/e3p from the exclusion plots obtained in the Bugey reactor experiment |29] and in 
the recent CHOOZ experiment [H3], which is the first reactor z/g — > i^e LBL experiment. 



At any fixed value of Arrig^^ in the range explored by the CHOOZ experiment we obtain 



the upper bound B, 
we have 



ee < B^^^ for a 



e, /i. From Eq.( [I.24j ), for the mixing parameter |[/, 



e3| 



\Ue,?<al 



or 



\U, 



e3| 



>l-a 



with 



1 - Bt 



;i.26) 



In Fig. |1.2| , we have plotted the values of the parameter ag obtained from the 90% CL 
exclusion plots of the Bugey and CHOOZ experiments. One can see that Cg is very small 
for Arnli > 10^^ eV^. Thus, the results of the reactor oscillation experiments imply that 
|f/g3p can only be small or large (close to one). 



Now let us take into account the results of solar neutrino experiments. The probability 
of solar neutrinos to survive in the case of a neutrino mass hierarchy is given by 



P::^.SE) = (1 - \Ue3\'Y P!,lXiE) + \Ue3\\ (1.27) 

where E is the neutrino energy and PuJ-JuX^^) is the two-generation survival probability 
of solar z/e's. If l^esP > 1 — ^e, from ( |1.27|) it follows that at all solar neutrino energies 
Pt^^u ^ 0.92. This is not compatible with the results of solar neutrino experiments. 
Thus, the mixing parameter |t/e3p must be small: 

\Uez?<al. (1.28) 

1.3.1 Decoupling of solar and atmospheric neutrino oscillations 

Let us consider the solar and atmospheric neutrino anomalies assuming that Arrigi is re- 
sponsible for the oscillations of solar neutrinos and Am|^ is responsible for the oscillations 
of atmospheric neutrinos (here we follow the discussion presented in Ref. ||46[|). 



From Fig. |L^ one can see that, if Amgj^ > 10~^ eV^ as indicated by the solution of the 
Kamiokande atmospheric neutrino anomaly, [|l^, the results of the CHOOZ experiment 
implies that 

If^esT < 5 X 10~^ (1.29) 

The averaged survival probability of solar electron neutrinos is given by Eq.( p..27| ). 
Taking into account the small upper bound (|1.29|) for |t/e3p, we have 

p::t.sE) ^ pll\{E) , (1.30) 

where PvJ-Jv^{E) is the two-generation survival probability of solar z/g's which depends on 

,2 „„i „. o _ \Ue2\ 



Arrig^jj = Am2i and sin^sun = ^^ , = - \Ue2\ ■ (1-31) 

43 



a-lf/,^12 



Hence the two-generation analyses of the solar neutrino data are appropriate in the three- 
neutrino scheme with a mass hierarchy and they give information on the values of 

Am2i = Am^^j^ and \Ue2\ — sin-i^sun- (1-32) 

The evolution equation for the flavor amplitudes ip^ (a = e, /x, r) of atmospheric 
neutrinos propagating in the interior of the Earth can be written as (see [^7| , |48[| ) 

t^^^ = ^{UM^U^ + A)m, (1.33) 

with 

^ = Vm , M^ = diag(0, Am^,, Am^,) , A = diag(Acc, 0, 0) , (1.34) 

and Ace = 2EVcc, where Vcc = V^GpNe is the charged- current effective potential 
which depends on the electron number density Ne of the medium {Gp is the Fermi 



constant and for anti-neutrinos Ace must be replaced by Ace = ~^cc)- If the squared- 
mass difference Am^i is relevant for the explanation of the solar neutrino problem, we 
have 



2B «1. (1.35) 

where i?© = 6371 Km is the radius of the Earth. Notice, however, that caution is needed 
for low-energy atmospheric neutrinos if Am^i ^ 10~^eV^, as in the case of the large 
mixing angle MSW solution of the solar neutrino problem and marginally in the case of the 
small mixing angle MSW solution. Indeed, ii Am^^ > 10~^ eV^ we have Am2iR^/2E <C 1 
only for E ^ 150 MeV. In this case, in order to get information on the three-neutrino 
mixing matrix with a two-generation analysis it is necessary to analyze the atmospheric 
neutrino data with a cut in energy such that Am2i-R©/2£^ -C 1. In order to be on the safe 
side, when the case of the MSW solutions of the solar neutrino problem are considered 
one can take into account the information obtained from the two-generation fit of the 
SuperKamiokande multi-GeV data alone |^ . 



The inequalities ( 1.35| ) imply that the phase generated by Am2i can be neglected for 



atmospheric neutrinos and M can be approximated with 

M^ ~ diag(0, 0, Am^^) . (1.36) 

In this case (taking into account that the phases of the matrix elements Ua3 can be 
included in the charged lepton fields) we have 

{UM^ U%>a ^ Ami, \Ua'3\ \Ua3\ • (1-37) 



Comparing this expression with Eqs.( |1.30D and (|1.32|) , one can see that the oscillations 



of solar and atmospheric neutrinos depend on different and independent Am^'s and on 
different and independent elements of the mixing matrix, i.e. they are decoupled. Strictly 
speaking \Ue2\ in Eqs.( |1.30D and ( |1.32| ) is not independent from \Ue3\ because of the 



unitarity constraint |t/eiP + I^e2p + l^^sP = 1? but the limit (|1.28|) on |?7e3p implies that 
its contribution to the unitarity constraint is negligible. 

Hence, we have shown that the smallness of |f/e3p inferred from the results of the 
CHOOZ experiment imply that the oscillations of solar and atmospheric neutrinos are 



decoupled [46 



From Eqs. (|1.33|) and (|1.37|) one can see that unless \Ue3\ -C 1, the evolution equations 



of the atmospheric electron neutrino amplitude ipe and those of the muon and tau neutrino 
amplitudes ipfi and ipr are coupled. In this case matter effects can contribute to the 
dominant u^ -^ Uj. oscillations (see [Q) and the atmospheric neutrino data must be 
analyzed with the three generation evolution equation (|1.33| ). 

From the results of the CHOOZ experiment it follows that the quantity |f/e3p is 
small and satisfy the inequality (|1.28|) (for Ami, ~ 10~^eV^). However, the upper 



bound for {Ue^l implied by Eq.( 1.28| ) is not very strong: If/gsl < 0.22. In the following 



we will assume that not only |f/e3p ^ 1? but also the element \Ue3\ that connects the 
first and third generations is small: |f/e3| *^ 1 (let us remind that in the quark sector 
2 X 10~^ < \Vub\ < 5 X 10"^). We will consider the other elements of the mixing matrix 
as free parameters and we will see that these parameters can be determined by the two- 
neutrino analyses of the solar and atmospheric neutrino data. In Section 1.3.3| it will 



10 



be shown that the hypothesis |?7e3| *^ 1 can be tested in future long-basehne neutrino 
oscillation experiments. 

If {Uesl ^ 1) for the evolution operator in Eq.( p..33| ) we have the approximate expres- 
sion 

/^ \ 

t/M^f/t + A^Am^J \U,s\' K-sWUrsl ] (1-38) 

V \Ur3\\U,,s\ \Ur3? J 

which shows that the evolution of z/g is decoupled from the evolution of u^ and z/^. Thus, 
the survival probability of atmospheric z/g's is equal to one and u^ -^ Vj. transitions are 
independent from matter effects and are described by a two-generation formalism. In this 
case, the two-generation analyses of the atmospheric neutrino data in terms of z/^ -^ v^ 
are appropriate in the three-neutrino scheme under consideration and yield information 
on the values of the parameters 

Amgi = Am^tn, and | f/^3 1 = sin t^atm ■ (1-39) 

1.3.2 The mixing matrix 

Under the assumption ^ei\ "^ 1, the values of all the elements of the three-neutrino 
mixing matrix can be obtained from the results of the two-generation fits of the data of 
solar and atmospheric neutrino experiments. The simplest way to do it is to start from 
the Chau and Keung parameterization of a 3 x 3 mixing matrix []TU|,|S^: 



C12C13 S12C13 Si3e"^i3 

U = \ -S12C23 - Ci2S23Si3e*'^i2 C12C23 - Si2S23Sl3e"'i=^ S23C13 | , (1.40) 

S12S23 - Cl2C23Sl3e*'^" -C12S23 - Si2C23Sl3e*'^i3 C23C13 

where Cij = cos ^ij and s^ = sin i)ij and 613 is the Dirac CP- violating phase (we have 
omitted three possible additional CP-violating phases in the case of Majorana neutrinos, 
on which there is no information). 

A very small |f/e3| implies that |si3| <^ 1. Since the CP-violating phase ^13 is associated 
with Si3, it follows that CP violation is negligible in the lepton sector^ and we have 



U ^ -S12C23 C12C23 S23 • (1.41) 



C12 


S12 


<] 


-S12C23 


C12C23 


S23 


■Sl2'S23 


-C12S23 


C23 



Using the information on |si2| — \Ue2\ and |s23| — \U,_i3\ given by the two-generation 
analyses of the results of solar and atmospheric neutrino experiments, for the moduli of 
the elements of the mixing matrix we obtain ||46|| : 



~1 0.03- 


-0.05 


<1 


0.02-0.05 0.71- 


-0.87 


0.49-0.71 


0.01-0.04 0.48- 


-0.71 


0.71-0.87 



Small mixing MSW: 0.02-0.05 0.71-0.87 0.49-0.71 , (1.42) 



^This can also be seen by noticing that the Jarlskog rephasing-invariant parameter |5l|~|53|] vanishes 
if one of the elements of the mixing matrix is zero. 
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Large mixing MSW: 



Vacuum oscillations: 



0.87- 


-0.94 


0.35- 


-0.49 


«1 \ 




0.25- 


-0.43 


0.61- 


-0.82 


0.49-0.71 


, (1.43) 


0.17- 


-0.35 


0.42- 


-0.66 


0.71-0.87 / 




0.71- 


-0.88 


0.48- 


-0.71 


«1 \ 




0.34- 


-0.61 


0.50- 


-0.76 


0.51-0.71 


• (1-44) 


0.24- 


-0.50 


0.36- 


-0.62 


0.71-0.86 / 





Let us remark that in the case of the small mixing angle MSW solution of the solar 
neutrino problem | f/gs I -C 1 could be of the same order of magnitude as | Ue2 \ ■ 

The mixing matrix ( |1.44| ) valid in the case of vacuum oscillations of solar neutrinos 



includes the popular bi-maximal mixing scenario that has been assumed in Refs. |54 



It is interesting to notice that, because of the large mixing of z/^ and z/,- with z/2, the 
transitions of solar z/g's in z/^'s and z/^-'s are of comparable magnitude. However, this 
phenomenon and the values of the entries in the (z/^, z/t-)-(z/i, z/2) sector of the mixing 
matrix cannot be checked with solar neutrino experiments because the low-energy z/^'s 
and z^T-'s coming from the sun can be detected only with neutral-current interactions, 
which are flavor-blind. Moreover, it will be very difficult to check the values of \U^i\, 
\Un2\, \Uti\ and \Ur2\ in laboratory experiments because of the smallness of 7712- 



In the derivation of Eqs.( |1.42D -( p..44| ) we have assumed that \Ue2\ < |f^ei| and |f/^ 



/x3 



< 



\Ur3\. The other possibilities, \Ue2\ > \Uei\ and \U^3\ > \Ur3\, are equivalent, respectively, 
to an exchange of the first and second columns and to an exchange of the second and 
third rows in the matrices (|1.42|) - (|1.44|) . Unfortunately, these alternatives are hard to 



distinguish experimentally because of the above mentioned difficulty to measure directly 
the values of \U^i\, \U^2\, \Uti\ and \Ur2\- Only the choice \Ue2\ < |t^ei|, which is necessary 
for the MSW solutions of the solar neutrino problem, could be confirmed by the results 
of the new generation of solar neutrino experiments (SuperKamiokande, SNO, ICARUS, 
Borexino, GNO and others ||55|) if they will allow to exclude the vacuum oscillation 
solution. 



1.3.3 Accelerator long-baseline experiments 

Future results from reactor long-baseline neutrino oscillation experiments (CHOOZ 
Palo Verde |S^, Kam-Land |^) could allow to improve the upper bound ( |1.28| ) on |f/e3p. 
In this section we discuss how an improvement of this upper bound could be obtained 
by future accelerator long-baseline neutrino oscillation experiments that are sensitive to 
^fi -^ ^e transitions (K2K [^| 



MINOS g, ICARUS g and others |TH6 



If matter effects are not important, in the scheme under consideration the parameter 



sin 2'due measured in z/,, 



'/xe 



long-baseline experiments is given by (see 



sm 



2-0 



fie 



4it/e3rit/^3r 



;i.45) 



If accelerator long-baseline neutrino oscillation experiments will not observe z/^ — > u^ 
transitions and will place an upper bound sin^ 2'i9^e < sin^ 2'&j^^^' , it will be possible to 
obtain the limit 



sin 2 27?^™''''^ 

l^^^l < nrr 12 

^I^M3l(min) 



:i.46) 
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where |C^/i3|Lin) is the minimum value of {U^^l"^ allowed by the solution of the atmospheric 
neutrino anomaly and by the observation of i/^ -^ Vr long-baseline transitions. For 
example, taking \U^z\l^^^ = 0.25 (see Eqs.(|Li2D-(|0§) we have |?7e3p < sin^ 2^^''^ 

If a value of sin^ 2'di^ ^ 10~^, that corresponds to the sensitivity of the ICARUS 
experiment for one year of running [Q, will be reached, it will be possible to put the 
upper bound jf/esl ^ 3 x 10~^. 

The observation of u^^ —>■ v^ transitions in long-baseline experiments will allow to 
establish a lower bound for |f/^3p because the parameter sin^ 2^^,- is given in the scheme 
under consideration by (see [^ ,|6^ 



sin^2^^. = 4|f/^3r|f/.3r- (1-47) 

From the unitarity relation |f/e3p + |f/^3p + |f^r3p = 1 it follows that an experimental 
lower bound siv?2§^r > siv?2'di^ allows to constraint |f/^3p in the range 

^1 - v/l-sin2 24r^^ < |f/,3p <\{^ + V^l-sin2 24r)) . (1.48) 

If sin^ 2'd^^T is found to be close to one, as suggested by the solution of the atmospheric 

neutrino problem, the lower bound |^^3|(inin) = \\^ ~ y ^ ~ sin^ 2'di^ J is close to 1/2. 

If matter effects are important, the extraction of an upper bound for |f/e3p from the 
data of z/^ — > z/g accelerator long-baseline experiments is more complicated. In this case 



the probability of z/^ — > z/g oscillations is given by (see [H3]) 



Ami-, I I eoi ^jni. 



'31 



(1.49) 

where E is the neutrino energy and L is the distance of propagation. This probability 
depends on the neutrino energy not only through the explicit E in the denominator of 
the phase, but also through the energy dependence of Aqc = '^EVqc- For long-baseline 
neutrino beams propagating in the mantle of the Earth the charged-current effective 
potential Vcc = V^GpNe is practically constant: N^ ~ 2iVyicm~^ {Na is the Avogadro 
number) and Vcc — 1-5 x 10^^'^ eV. 

If long-baseline experiments will not observe z/^ -^ Ve transitions (or will find that 
they have an extremely small probability) for neutrino energies such that Ace ^ ^^^si; 
it will mean that |f/e3p is small and a fit of the experimental data with the formula 
(|1.49|) will yield a stringent upper limit for |t/e3p (taking into account the lower limit 



|f^/^3p > I^Mslfmin) obtained from the solution of the atmospheric neutrino anomaly and 
from the observation of v^ -^ v^ long-baseline transitions). On the other hand, the non- 
observation of z/^ — > z/g transitions for neutrino energies such that Ace ^ ^'^si does not 
provide any information on |t/e3p because in this case the transition probability (|1.49| ) is 
suppressed by the matter effect. Hence, in order to check the hypothesis \Uez\ <^ 1, as 
well as to have some possibility to observe z/^ — *> z^e transitions if this hypothesis is wrong, 
it is necessary that a substantial part of the energy spectrum of the neutrino beam lies 
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below 



Amii 



2V( 



30GeV 



cc 



10-2 eV^ 



:i.50) 



This requirement will be satisfied in the accelerator long-baseline experiments under 
preparation (K2K |58|, MINOS [||], ICARUS ^ and others |l|-|l) if Amii is not 
much smaller than 10~2eV^. 



1.3.4 Neutrinoless double-/? decay 



Let us consider now the implications of the result in Eq. ( |1.28| ) for neutrinoless double-/? 
decay experiments. The investigation of neutrino oscillations does not allow (see P,|^ 
to answer the fundamental question: are massive neutrinos Dirac or Majorana particles? 
Only investigations of neutrinoless double-/? decay could allow to answer this question. In 
the case of a three-neutrino mass hierarchy for the effective Majorana mass we have [S3,BB[] 



|(m)|^|f/,3p/ 



Amii 



;i.51) 



Taking into account the bound (|1.28|) on |t/e3p; we obtain the following constraint for 
the effective Majorana mass in neutrinoless double-/? decay 



IWI 



< 



Amii 



(1.52) 



The value of this upper bound as a function ATTig^ obtained from 90% CL exclusion 



plots of the Bugey [^ and CHOOZ ^^ experiments is presented in Fig. p..3| (the solid 
and dashed line, respectively). The region on the right of the thick straight solid line is 
forbidden by the unitarity bound | (m) | < \J l\m^^. 

Also the results of the Super-Kamiokande atmospheric neutrino experiment |T^ imply 
an upper bound for |f/e3p. The shadowed region in Fig. p..3| shows the region allowed by 



Super-Kamiokande results at 90% CL that we have obtained using the results of three- 
neutrino analysis performed by Yasuda |6^ . 

Figure |1.3| shows that the results of the Super-Kamiokande and CHOOZ experiments 
imply that |(m)| < lO^^eV. 

The observation of neutrinoless double-/? decay with a probability that corresponds 
to a value of |(?7t.)| significantly larger than 10~2eV would mean that the masses of three 
neutrinos do not have a hierarchical pattern and/or exotic mechanisms (right-handed 
currents, supersymmetry with violation of R-parity, . . . |^,^) are responsible for the 
process. 

Let us notice that from the results of the Heidelberg-Moscow ^^Ge experiment |70| 
it follows that \{v(i)\ < 0.5 — 1.5 eV. The next generation | 
experiments will reach \{vti)\ ~ 10~^eV. Possibilities to reach |(m)| 



of neutrinoless double-/? 
10~2 eV are under 



discussion |^| 



1.4 Discussion 

Xing: A hierarchy of Am^'g does not imply a hierarchy of masses. 
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Bilenky: Yes, I agree. However, on the basis of a similarity with the observed hierarchies 
of the masses of quarks and charged leptons, a hierarchy of Attic's is a strong hint 
in favor of a hierarchy of neutrino masses. 
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2 A review of recent neutrino oscillation solutions 
and their implications for the future experimental 
neutrino program — Achim Geiser 

2.1 Introduction 

Currently there is a wealth of unexplained phenomena in neutrino physics which can be 
interpreted as indications for the existence of neutrino oscillations. The solar neutrino 
problem [jl|, lacking a satisfactory astrophysical solution, is commonly attributed to the 
disappearance of electron neutrinos into some other neutrino type. In the wake of the 
recent Super-Kamiokande results [0] , the atmospheric neutrino anomaly [Q] is interpreted 
as evidence for neutrino oscillations involving z/^ disappearance. The LSND experiment [^ 
claims direct evidence for u^ — z/g oscillations in a region which is partially unconstrained 
by other experiments. Furthermore, neutrinos are a prime candidate for a partial solution 
to the missing dark matter problem [^ if at least one mass eigenstate lies in the eV range. 
A more detailed discussion of the experimental indications will be presented elsewhere. 

It is the purpose of this contribution to systematically review the different classes of 
possible solutions and discuss some of their features, including their predictions for future 
experiments. An attempt is made to include all oscillation scenarios published after 1995, 
and relevant preprints not older than one year. First, all these scenarios will be classified 
according to general criteria. Then, some of these models will be discussed in more detail. 
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2.2 Classification of neutrino oscillation solutions 

Results of neutrino oscillation experiments are often expressed in terms of an effective two 
flavour oscillation scheme with a mixing angle sin^ 26 between the two flavours and a mass 
difference 5m? between the two relevant mass eigenstates. The indications for neutrino 
oscillations from solar neutrinos, atmospheric neutrinos and LSND each suggest a different 
value of dm?. Unfortunately it is not possible to find a unique solution simultaneously 
satisfying all these indications. 

Solutions involving the 3 known active neutrino flavours yield only two independent 
5m? values. The problem is therefore overconstrained, and some experimental evidence 
has to be discarded (or equivalently, a bad fit has to be accepted) in order to find a 
solution. This induces some arbitrariness in which parts of the data are deemed to be 
reliable, and which parts should be ignored. Once some information has been discarded, 
the remaining information is often not sufficient to uniquely constrain the 3 mixing angles. 

Solutions involving 1 or more sterile neutrinos are generally underconstrained due to 
the 6 or more mixing angles involved. Therefore many of the models of this kind make 
some simplifying assumptions which reduce the parameter space before attempting to find 
a solution, or impose constraints obtained from specific GUT models or other extensions 
of the standard model. 

In general, there are four big classes of possible solutions 

• Solutions trying to accommodate all the experimental evidence in the standard 3 
neutrino scheme, at the expense of accepting a bad fit to part of the data |p|-P . 



Solutions which simply discard one of the experimental indications p^-p3[|. A good 
fit to the remaining indications can then be obtained from three active neutrinos. 
The current "standard" solution (dropping LSND) falls into this class. 

Solutions which invoke 1 additional light sterile neutrino. This is the minimum 
required to accommodate the three different 5m? values discussed above. Commonly 
the sterile neutrino is used to solve the solar neutrino problem p4|-p^ , but scenarios 



explaining the atmospheric neutrino anomaly through active-sterile oscillations are 
equally possible p7|^9|. Alternatively, the sterile neutrino can be used to explain 
LSND |^0|- Many models predicting such scenarios have been investigated |3l[| . 



Solutions which invoke more than one sterile neutrino (usually 3, motivated by 
assuming some symmetry between the active and sterile neutrinos) ||3^36|. It has 



recently been shown that Big Bang nucleosynthesis limits which seemed to exclude 
some of these scenarios can be evaded 



Other solutions can be derived by partially replacing the oscillation hypothesis by 
other non-standard model effects (neutrino magnetic moments, anomalous interactions 
...) or by discarding more of the experimental evidence than the minimum needed 



to satisfy a given scenario. 

In order to put some structure into the discussion of the multitude of oscillation 
scenarios, we adopt a formal classification of each solution according to its interpretation 
of the various pieces of experimental evidence. 

The classification criteria used are 
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A. The solar neutrino problem is 

• I. not caused by neutrino oscillations (astrophysical solution, spin flip due to neu- 
trino magnetic moment, ...) 

• II. due to z/g -^ ^fi (or u^) oscillations 

• III. due to z/g -^ Us oscillations 

B. The atmospheric neutrino anomaly is 

• I. not caused by neutrino oscillations (anomalous Ur interactions, ...) 

• II. due to u^ -^ Ur oscillations 

• III. due to a linear combination of v^ — Ve and v^ — »• Vr oscillations 

• IV. due to z/^ -^ vs oscillations 

C. The LSND result is 

• I. not caused by neutrino oscillations 

• II. due to direct v^ -^ v^ oscillations (effective 2x2 mixing matrix) 

• III. due to indirect z/^ -^ z/g oscillations (full 3x3 matrix is relevant) 

In addition, all scenarios can be classifled in terms of the different MSW and vacuum 
oscillation schemes used to solve the solar neutrino problem [Q. Further classification 
criteria could be the choice of the relevant (5m^, the level of compatibility with existing 
reactor and accelerator limits and limits on double beta decay, the compatibility with 
different dark matter scenarios, and the compatibility with Big Bang nucleosynthesis and 
supernova processes. An explicit formal consideration of these criteria would however be 
impractical to handle. They will therefore be discussed elsewhere. 

Table |2.1| classifies all recent neutrino oscillation solutions according to categories A.- 
C, and lists the corresponding expectations for current and future neutrino oscillation 
experiments. It is interesting to note that there is almost no combination of elements 
of this classification scheme which has NOT been discussed in the recent literature as a 
possible scenario. Conversely, there are only very few linear combinations of results of 
future measurements which would NOT correspond to a possible oscillation scenario. 

2.3 Discussion of selected models 

2.3.1 Three neutrino solutions including LSND 



Before the recent results from CHOOZ [Q and Super-Kamiokande [0], this used to be one 
of the favorite options. Discarding the solar neutrino deficit, one can easily build models 
including LSND and atmospheric neutrinos only |T^,^. However, currently no serious 



alternative to the neutrino oscillation hypothesis of the solar neutrino deficit exists. On 
the other hand, one could discard the atmospheric neutrino anomaly and build models 
using LSND and solar neutrinos only [|1^ . Going one step less far, the model of Cardall 
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Table 2.1: Classification of neutrino oscillation solutions in terms of their interpretation of the current 
experimental data. The predictions for selected ongoing and future experiments are given in each case. 
The roman numbers refer to the classification scheme described in the text. The symbol + stands for 
the expectation of a positive effect (deviation from the no oscillation expectation) and implies that the 
corresponding scenario is excluded if no signal is observed. The symbol - indicates the expectation of 
a negative result and implies that the scenario is excluded if a positive signal is seen. The symbol 
indicates that both positive and negative results are possible in a given scenario, often depending on the 
results of other future measurements. In many cases it means that the scenario is favoured in the case 
of a po sitiv e effect, but not (completely) excluded if no signal is seen. The experiments are explained in 
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and Fuller tries to reconcile LSND, atmospheric, and solar neutrinos by ignoring the 



energy dependence of the up-down asymmetry in the Kamiokande Multi-GeV data p4 



The same dim? can then be used for LSND and atmospheric neutrinos. This scenario 
is now strongly disfavored by the much more significant up/down asymmetry observed 
by Super-Kamiokande p[. An alternative is to merge the 5w?' of the atmospheric and 



solar neutrinos by discarding the Homestake |45| result or by somewhat increasing its 
error |iB|. This allows an energy- independent suppression of the solar neutrinos at high 
Srn?. However, the corresponding model of Acker and Pakvasa ||^ is not consistent 
with the recent CHOOZ and Super-Kamiokande results, which strongly constrain the 
possible u^ — z/g contribution [^. A new model by Thun and McKee merges the 
Cardall/Fuller and Acker/Pakvasa schemes by letting both the solar neutrino and the 
LSND 5im? contribute significantly to the atmospheric neutrino anomaly. This model 
can successfully describe the Super-Kamiokande up/down asymmetry of the e//i ratio, 
but fails to fit the shape of the individual v^ and Ve fluxes and asymmetries. It is therefore 
excluded unless large unknown systematic effects are assumed (i.e. the corresponding flux 
measurements systematically disagree ||5^). On the other hand, this kind of model might 
be interesting in the context of a recent claim that the KARMEN anomaly |^8[ can be 
explained through neutrino oscillations [HO]. 



2.3.2 Three neutrino solutions excluding LSND 

One currently favoured way out of the dilemma is to discard the LSND result 0], which 



has so far not been confirmed by KARMEN |^ or other experiments [^ . This leads to 



a scenario which is very close to one of the standard scenarios before the LSND claim |T2 



The atmospheric neutrino anomaly is interpreted in terms of v^ — v^ oscillations, while the 
solar neutrino deficit is ascribed to v^ — u^ oscillations. Both MSW and vacuum solutions 
are allowed in this context |l[]. Variations of this scheme allow an additional significant 
^iJL~i-'e contribution to the atmospheric neutrino anomaly |T^,|n[. In particular, the three- 
fold maximal mixing scenario |jl6[ often cited as model of Harrison, Perkins, and Scott 
is not completely excluded by CHOOZ and Super-Kamiokande [O. However, the pure 



— Vf. interpretation of the atmospheric neutrino data is no longer allowed |2^. Models 



explaining both the solar and atmospheric neutrino problems through i/^ — z/g oscilla- 
tions [|2T| are therefore excluded. Motivations why neutrino mixing angles could/should 



be large are given by various authors [^^^H • 



2.4 Sterile neutrinos 

If the current indications for three independent dirP' are confirmed, the only way out 
is the introduction of at least a fourth neutrino. LEP [^] has measured the number 
of neutrinos coupling to weak interactions to be N^^ = 2.994 ± 0.012, corresponding to 
the three known lepton generations. Any extra neutrinos must therefore be either very 
massive (rriu > Mz/2) or sterile with respect to weak interactions. In a minimal extension 
of the standard model, these could e.g. be the right-handed (left-handed) partners of the 
standard model neutrino (antineutrino). A tiny lepton-number violating interaction could 
then induce neutrino-antineutrino oscillations, in analogy to K^ — K^ oscillations |53 



Since oscillations do not change the spin orientation, left-handed neutrinos would oscillate 
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into left-handed antineutrinos which appear sterile. Alternatively, extra neutrino singlets 



or multiplets can be introduced in the context of Grand Unified Theories (see e.g. p9| 



and references therein) or other extensions of the standard model |^3|]. Many of these 



schemes invoke the see-saw mechanism |Q , which makes some of the additional neutrinos 



too heavy to be detected. Current experiments can not clearly distinguish oscillations 
with sterile neutrinos from standard flavour oscillations for either atmospheric [^] or 
solar 1^ neutrinos, u — Ug oscillations are therefore allowed in both cases. Finally, it 
is worth noting that many of the models advocating light sterile neutrinos are justified 
without the explicit requirement of the LSND constraint, and hence do not depend on 
LSND being confirmed. 

2.4.1 Four neutrino solutions 

Many models consider one sterile neutrino in addition to the active ones. This is the 
minimum required to obtain three independent 5m^'s, and is sometimes motivated by 
specific GUT models. Two classes of such models can be distinguished, depending on 
which of the experimental observations is ascribed to active/sterile oscillations. Schemes 
invoking z/g — Ug oscillations for the solar neutrino deficit ||2^,|25| are the oldest class of 
such models. In this context, the atmospheric neutrino deficit is explained by z/^ — u^ 
oscillations, and LSND by u^ — u^, oscillations. More complicated variants of these models 



are possible p6|,^. Alternatively, v^ — Pg oscillations can be invoked to describe the 
atmospheric neutrino anomaly p7|-p9|. In this case, the solar neutrino deficit can be some 
linear combination of v^ — i^^i, ^e — ^^t, and v^ — Vs oscillations. The possibility that the 



LSND signal is due to indirect oscillations involving the Vg has also been considered pO 



Mass textures for all these solutions are being investigated |31 



2.4.2 Five or six neutrino solutions 

Models with more than four neutrinos have been unpopular due to their apparent conflict 



with bounds on the number of neutrinos from Big Bang nucleosynthesis p5| , ^ . However, 
it has recently been demonstrated that these limits can be evaded through modifications 
of the Big Bang nucleosynthesis model which are a consequence of the active/sterile 



neutrino oscillations themselves |^. Scenarios with five light neutrinos are considered 



in [^. A more natural configuration of 6 neutrinos can be obtained by assigning a light 
sterile (right-handed) neutrino to each lepton generation through some symmetry (e.g. 
parity). In the context of such scenarios |3^,|3^, both atmospheric and solar neutrino 



oscillations are explained through close to maximal active/sterile oscillations, while the 
LSND result is interpreted as standard v^ — v^ oscillations. Other phenomenologically 
similar models |^| invoke an unstable z/^ in the MeV mass range to evade the Big Bang 



nucleosynthesis constraint. 



2.5 Implications for future experiments 



As can be inferred from table ^3| almost any linear combination of the outcome of ongoing 
and future experiments corresponds to a possible oscillation scenario. The predictions 
for these experiments therefore depend on the experimental and theoretical bias of which 
current experimental results are believed to be true and how they are interpreted. 
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No single experiment can uniquely conGrm any of the discussed scenarios, but certain 
combinations of future experimental results will rule out or support some of them. (Mini)- 
Boone |5^ and/or PS 1216 |]58|, if approved, will directly check the LSND claim, therefore 
ruling out half of the listed scenarios. Similarly, the measurement of the NC/CC ratio 
at the Sudbury Neutrino Observatory (SNO |5^) will uniquely check the z/e — z/^ (or 
z/g — Ut-) interpretation of the solar neutrino deficit. Measurements of the NC/CC ratio 
in atmospheric neutrinos, for instance by Super-Kamiokande [^] are very powerful in 
principle for the distinction of atmospheric u^ — z/^/g and z/^ — Ug oscillations. In practice, 
they might be limited by low statistics and/or large systematic errors |^T|. Alternatively, 
the combination of future long [|62|,|63| and short [^ or intermediate ||6 



baseline r 
can achieve 
the same distinction with much higher significance, provided the atmospheric 6m'^ is not 



appearance experiments (but not the long baseline experiments alone |61 



too small. With the same caveat, long baseline disappearance experiments []63|, |66|, |67|1 
will check the atmospheric neutrino anomaly, and appearance experiments p2 



3 will 

test any significant z/^ — z/g contribution. The observation of an oscillation signal beyond 
the current limits at the CHOOZ [^] or Palo Verde ||B^ reactors would establish such 
a contribution. The observation of u^ — z/,- oscillations in ongoing (CHORUS/NOMAD 

short or intermediate baseline experiments could force the z/^ — Ug 



or future 



interpretation of the atmospheric neutrino anomaly. Alternatively, it could confirm the 



oscillation scenario II for LSND ||70[ invoked in the Cardall/Fuller, Acker/Pakvasa, and 
Thun/McKee schemes |]7T|. The KamLand experiment ||72[| will check the large angle 
MSW solution of the solar neutrino problem |jl| , while the Iodine [^ and Borexino |7^ 



experiments will try to confirm the energy dependence of the solar neutrino deficit. 

2.6 Conclusions 

An attempt has been made to classify all recently proposed neutrino oscillation solutions 
according to their interpretation of the current experimental evidence, and to discuss 
their implications for future experiments. At present, no single scenario is emerging 
as "the" obvious candidate. Essentially all planned future neutrino experiments are 
needed and useful to constrain and/or exclude certain classes of scenarios, and almost any 
combination of future experimental results corresponds to a possible oscillation scheme. 



2.7 Discussion 

Fogli: A detailed quantitative calculation shows that the model of Thun and McKee 
(as well as the other models in reference ||^) gives a very bad fit of the sepa- 
rate zenith angle distributions of atmospheric e-like and /x-like events measured by 
Super-Kamiokande (see Figs BTT] and p.2[). The model of Thun and McKee assumes 



that the atmospheric neutrino anomaly is explained by dominant z/^ -^ z/g transi- 



tions with Arrigi 



< 



10-'^ eV' 



I.e. 



43 . Sub-dominant 



u. 



oscillations due to Am^a ~ 



below the CHOOZ bound 
0.4 eV^ cause an additional energy- independent 
suppression of the atmospheric muon neutrino flux. 

As one can see from Fig. |2.1| (obtained with the method described in the second pa- 
per of Ref. |T2[)) the resulting distribution of /i-like events in the Super-Kamiokande 
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detector do not fit the data because they are too fiat. The reasons for this flatness 
are: 

1. The matter effect suppresses the effective mixing. From the comparison of 
Figs |2.1| and \2.2\ one can see that the matter effect is not neghgible (a detailed 
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calculation shows that the matter effect can give up to 60% change of the /x/e 
ratio of upward-going events). 

2. The presence of the sub-dominant energy-independent u^ -^ Vj. oscillations. 

Furthermore, one must take into account that the model of Thun and McKee gives 
a bad fit of solar neutrino data and a marginal fit of reactor and accelerator neutrino 
oscillation data. 

Geiser: I verified that it gives a bad fit of the separate e-like and /i-like Super- 
Kamiokande data as a function of L/ E. However, this incompatibility might be 
attenuated by changing the spectral slope of the cosmic rays used for the theoret- 
ical calculations of the atmospheric neutrino fluxes, (note added: after studying 



figs. |2.1j and ^^ in more detail, I agree that it seems indeed hard to account for the 



Multi-GeV muon up-down asymmetry once matter effects are taken into account.) 

Frekers: Which are the astrophysical constraints on the scheme with six pseudo-Dirac 
neutrinos? In particular, are there constraints from Big-Bang nucleosynthesis and 
supernovae? What about the energy loss in supernovae? 

Geiser: The six pseudo-Dirac neutrino scheme [|^ is indeed compatible with Big Bang 
nucleosynthesis through the mechanism proposed by Foot and Volkas (Ref. |^). I 
did not study the compatibility with supernova processes in detail, but I am not 
aware of any incompatibility with cases treated in the literature so far. 

Langacker: Maximal mixing oscillations of active into sterile neutrinos in the supernova 
core may cause a too large energy loss. 

Giunti: The maximal mixing oscillations in the supernova core of active into sterile 
neutrinos is suppressed if the effective potential is bigger than all relevant Am'^. 
In this case the effective mixing angles are close to 7r/2, i.e. in practice there is 
no mixing. However, transitions of active into sterile neutrinos can occur while the 
neutrinos escape the supernova and should be studied in detail. 

Sarkar: Big Bang nucleosynthesis is well studied and certainly rules out 3 additional 
maximally mixed sterile neutrinos. With regard to evading this bound through 
the Foot- Volkas mechanism ^^ , please note that the sign of the lepton asymmetry 
generated is arbitrary so the value of Ni, can go up as well as down! 

Frekers: People proposing new models should explain clearly their compatibility with 
physical and astrophysical data. 
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Detection of z/^ -^ v^ at Am^ > IQ-^ eV^ in ICARUS 
— Stan Otwinowski 



The recent results of the Super-Kainiokande experiment appear to provide strong evidence 
for v^ —>■ Vx atmospheric neutrino oscillations [0]. The most interesting possibility is that 
the Vx is a z/,-. A large effort at FNAL and the CERN-Gran Sasso Laboratories will be 
directed at the detection of a r appearance. However if the relevant Am^ is less than 
10^'^ eV^ this is a very difficult task since the r threshold energy of 3.5 GeV forces the use 
of relatively high energy neutrinos. It is conceivable that the Super-Kamiokande results 
could be confirmed without determining the actual nature of the Vx (K2K or MINOS 
or ICARUS) leaving open the possibility that Vx is a sterile neutrino. There is another 
way to have some hint that Vx is an active neutrino by exploiting the possibility of three 
neutrino fiavour mixing. The basic idea is that at the same level if the Super-Kamiokande 
results are due to u^ -^ u^. then the processes z/e -^ u^ and u^ -^ z/^ may occur at the 
same Am"^. Observation of the z/^ -^ z/g process at foreseen electron energy region does 
not rule out the sterile neutrino but, we believe, the sterile neutrino hypothesis would be 
less favoured. 

We believe that ICARUS [0] is the ideal detector for this observation, since the excel- 
lent vertex resolution can be used to discriminate against neutral current n^ events. 

There are two possibilities for this detection: 

1. z/^ — i> z^e has a transition probability bigger than 10% (up to the CHOOZ limit |^). 
In this case the signal would be much larger than the u^. background in the LBL 
beam (0.8%) and detection is likely. It seems very possible that this will occur since 
the large z/^ — > z/^ transition probability suggests that the neutrino sector is unlike 
the quark mixing sector with the small CKM mixing angles — all of the mixing 
angles could be relatively large. 

2. The transition probability is less than 10^^ (say 5 x 10'^). It may still be possible 
to detect the z/^ -^ v^. transitions through the observation of a modification of the 
z/e energy spectrum. This is due to the fact that the normal z/g beam is somewhat 
harder than the v^ beam because of the z/g origin from K^ decays and the very small 
forward angles accepted by the detector in the LBL beam. For example, for certain 
values of the mixing parameters we expect considerable z/^ — > z/g transitions at 5 — 10 



GeV neutrino energy — at lower energy than the normal z/g beam (Fig. p.l| ). A 
lower energy v^ beam, as recently suggested by A. Ereditato in Ref. |^, would be 
even more useful (see Fig. pTl|). 
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As an example we consider a model of Cardall, Fuller and Cline 
analysis was studied. In the case of the mixing matrix 



U, 



CF,lma 
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where a 3-z/ 



(3.1) 



we have 



v„ 






p 
p 



1, 

10-1 



(3.2) 
(3.3) 



This model also fits the Solar neutrino data with the large mixing angle (Ima) solution 
at /\m? ^ 2 X 10~^ eV^. This is only an example and a more detailed analysis of the 3-z/ 
mixing is called for (see Ref. [Q). 



3.1 Discussion 

Langacker: I think that the observation of long-baseline z/„ — > Ue oscillations cannot 



help to distinguish the i/^ 
anomaly. 



Pt- and Pu -^ ^s solutions of the atmospheric neutrino 
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Giunti: I also think so. For example, with four neutrinos there can be simultaneous 
large i/^ -^ Ug oscillations and observable z/^ -^ u^ oscillations at the atmospheric 
scale of Am'^. However, in the case of four neutrinos large long-baseline z/^ — *> u^ 
oscillations are incompatible with the LSND indication in favor of short-baseline 
i^^l —^ ^e and Vy, -^ z/g oscillations 0. 
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4 Phenomenological Ansatze for Large Mixing of 
Three Neutrinos — Zhi-Zhong Xing 

4.1 Introduction 

Recently the Super-Kamiokande Collaboration has reported new and stronger evidence 
for the existence of the atmospheric neutrino anomaly. The data particularly favor an 
interpretation of the observed muon-neutrino deficit by z/^ <-^ Ur oscillations with the 
mass-squared difference Aml^^ Ri (0.5 — 6) x 10^^ eV^ and the mixing factor sin^ 2^atm > 
0.82 at the 90% confidence level |jl|]. The long-standing solar neutrino deficit has also 
been confirmed in the Super-Kamiokande experiment. Analyses of the energy shape 
and day-night spectra of solar neutrinos seem to favor the "Just-So" mechanism, i.e., 
the long- wavelength vacuum oscillation with Arn^^^ ^ 10^^° eV^ and sin^ 26'sun ~ 1 0- 
However, it remains too early to justify whether the solar neutrino anomaly is attributed 
to the "Just-So" oscillation or to the the matter-enhanced oscillation (i.e., the Mikheyev- 
Smirnov-Wolfenstein (MSW) mechanism). An analysis of all available solar neutrino 
data based on the latter gives two different ranges of the oscillation parameters [§ : the 
large-angle MSW solution with Am^^^ ^ 10"^ eV^ and sin^ 26'sun ~ 0.8 @] or the small- 
angle MSW solution with Am'^^^ ^ 5 x 10~^eV^ and sin^ 26'sun ~ 5 x 10"^. If the large 
mixing angles 6'atm and 6'sun in the "Just-So" or the large-angle MSW scenario are finally 
confirmed, one would have an indication that the physics responsible for neutrino masses 
and lepton fiavor mixing might be qualitatively different from that for the quark sector. 
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The small mixing angle ^sun appearing in the small-angle MSW scenario, on the other 
hand, seems similar to the small angles observed in the quark flavor mixing phenomenon. 

The LSND evidence for i/^ — >■ z/g oscillations [^], whose parameters lie in the ranges 
Attt-lsno ~ (0-4 — 2) eV^ and sin^ 26'lsnd ~ 10"^ — 10~^, was not conflrmed by the recent 
KARMEN experiment 0. Because a further examination of the LSND result will be 
available in the coming years, the conservative approach is to put it aside tentatively 
from the evidence of solar and atmospheric neutrino oscillations. Indeed it is extremely 
difficult, if not impossible, to accommodate all solar, atmospheric and LSND oscillation 
data within the scheme of three active neutrinos. 

It is remarkable that the CHOOZ experiment 0, in which the survival probability 
of Vf, neutrino was measured, indicates that sin^ 20ch < 0.18 if Am^y > 9 x 10""^ eV^. 
This result supports that from the Super-Kamiokande experiment, i.e., the atmospheric 
neutrino deficit comes most likely from v^ ^^ v^ oscillations instead of v^ <->■ t/g oscillations 
in the three-neutrino framework. In particular, it turns out that the (1,3) element of the 
3x3 lepton fiavor mixing matrix is naturally small in magnitude and the atmospheric 
and solar neutrino oscillations approximately decouple - they are separately dominated 
by a single mass scale (see, e.g., Refs. P,^): 



v, 



oscillation : Am?„„ = Amn, 



^ ^^ i^e ^aK.LLLa.hL\jLL . l^iiu^^^ — L^nu2i , 



T-'lj. *-^ ^T oscillation : Am^^^jj^ = Amgg ~ ^"^31 ) (4-1) 

where Am?- = |m^ — rp?A and mj (for % = 1,2,3) is the neutrino mass. The smallness 
of the (1,3) mixing element in this simple picture, compared with the smallness of the 
(1,3) element in the 3x3 quark mixing matrix [1^, refiects some kind of similarity 



between lepton and quark fiavor mixings. Whether our present understanding of solar 
and atmospheric neutrino oscillations as described in Eq. (1) is close to the truth or not, 
however, remains an open question. 

To interpret current data on atmospheric and solar neutrino oscillations, one may 
follow a phenomenological way to construct the lepton fiavor mixing matrix which can 
account for the experimentally favored values of 6'atm and ^sun- Any insight into the dy- 
namics of lepton mass generation, however, requires nontrivial theoretical steps beyond 
the standard model [|ll|. Note that the oscillation parameters Am? and sin^26' are in 
general expected to correlate with each other, as the latter comes from the mismatch 
between charged lepton and neutrino mass matrices and should depend on the ratio(s) 
of neutrino masses. One can therefore distinguish between two different types of models, 
in which sin^26' rely specifically on Am?, or in which sin^26' is independent of any neu- 
trino mass. In all experimental analyses, of course, Am"^ and sin^ 26 are treated as two 
independent parameters. 

Before a successful theory of lepton mass generation and fiavor mixing appears, the 
proper phenomenological approach might be first to identify the patterns of lepton mass 
matrices guided by some kinds of possible symmetries, and then to calculate the lepton 
fiavor mixing matrix. The "success" of such lepton mass and mixing ansatze can only be 
"justified" by their consequences, i.e., if they agree well with current data on neutrino 
oscillations and if they are consistent with some other constraints on neutrinos (e.g., the 
constraints from the neutrino less /3/3-decay or from the hot dark matter). We hope that 
our phenomenological attempts can provide useful hints towards the theoretical solution 
of lepton mass and mixing problems. 
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Following the strategy outlined above we are going to discuss a simple phenomeno- 
logical model of lepton masses and flavor mixing, on the basis of flavor democracy for 
charged leptons and mass degeneracy for neutrinos ]I2[ . We concentrate on three different 



symmetry-breaking scenarios in the scheme of this model, which can lead to the "nearly 
bi-maximal" mixing |T2|,|13|, "small versus large" mixing ||H| and "exactly bi- maximal" 



mixing |T5[ for solar and atmospheric neutrino oscillations, respectively. Implications of 
these ansatze, in particular on the upcoming long-baseline neutrino experiments, are also 
discussed. 



4.2 Three-neutrino mixing 



The hierarchy of ^m^i and Arrigg (or Ar/ig^) set in Eq. (1) has little implication on 
the relative magnitude of ttii, m2 and m^. For this reason, one has to assume speciflc 
neutrino mass spectra in constructing phenomenological models of lepton flavor mixing. 
Two distinct possibilities are of particular interest and have attracted a lot of attention: 

(a) Three neutrino masses perform a clear hierarchy: mi -C m2 *^ "^3 ; 

(b) Three neutrino masses are almost degenerate: nii k. 7712 ~ m^ . 

In comparison with the well-known mass spectrum of the charged leptons {rrie -C m^ -C 
rrir), we expect that the charged lepton and neutrino mass matrices {Mi and M^) in 
scenario (a) may both take the hierarchical form, similar to the quark mass matrices. 
The neutrino mass matrix in scenario (b), however, must take a form quite different 
from the charged lepton mass matrix. Therefore large mismatch between Mi and M^, in 
scenario (b), which leads in most cases to large flavor mixing, seems very natural. In 
addition, scenario (b) is welcome by the interpretation of hot dark matter, if the sum of 
three neutrino masses amounts to a few eV 
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It proves convenient to introduce a speciflc parameterization for the 3x3 lepton flavor 
mixing matrix V, which can be obtained from diagonalization of Mi and M^, in a chosen 
basis of flavor space. In view of our wealthy knowledge on quark mass matrices and quark 
flavor mixing, we expect that an appropriate description of lepton flavor mixing, in terms 
of three Euler angles 



V 




and one CP- violating phase 0, takes the form |13tp!7 










c 



-c„s 




(4.2) 



where s^ = sin 61, s^, = sin O^,, c = cos 6, etc. Throughout this work we assume CP invari- 
ance in the lepton sector, i.e., we take = 0. The mixing angles of this parameterization 
may have very instructive meanings, as one can see later on. In correspondence with 
the one dominant mass scale approximation made in Eq. (1) for solar and atmospheric 
neutrino oscillations, |^;| -C 1 is naturally expected. Then the mixing factors sin^ 2^atm 
and sin^ 29^un turn out to be 



sin 26'atm 

sin^ 26's„r, 



sin^ 26 , 
sin^ 26^ 



(4.3) 
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to a good degree of accuracy. 

It is not difficult to understand the smallness of the mixing angles 6i, which conies 
primarily from mixing between the ffist and second families of the charged lepton mass 
matrix M/. In some reasonable Fritzsch-like ansatze of lepton mass matrices (with or 



without the see-saw mechanism [0,|T9[), 6i is essentially given by 



Oi = arctan ( . /— j ^ 4° . (4.4) 

Therefore the decoupling condition for the solar and atmospheric neutrino oscillations 
can be satisfied. It should be noted, however, that non-vanishing 6i may have significant 
effect on transitions of z/g and Ue neutrinos in the long-baseline neutrino experiments. 

For example, one may carry out a long-baseline (LB) experiment in which neutrino 
oscillations are governed only by the mass scale Am^tm = ^7/132 ~ Ar7i|^ ^ Arn^i- In 
this case the transition probability of z/^ to Ue reads 

P{i^^^i^e)LB = sm''29ism^9 sin^ ('l.27^^^'] , (4.5) 

where P denotes the momentum of the neutrino beam (in unit of GeV), and L is the 
distance between the neutrino production and detection points (in unit of km). We 
observe that sin^2^;, like sin^ 26' and sin^ 2^j, in Eq. (3), does get an apparent meaning 
from neutrino oscillations. The transition probability of u^ to z/,- in the similar long- 
baseline experiment is given by 

Pif^e -^ '^r)LB = (l + taU^ Oi) COt^ 6 P{u^ ^ U,)^B ■ (4.6) 

As a result, measurable signals in Eqs. (5) and (6) (or one of them) in the future 
experiments will definitely rule out the "exactly bi-maximal" mixing ansatz of three 
neutrinos [|T^], in which 6i = holds. 

The question is now of whether the mixing angle 6^, depends on the neutrino mass 
ratios. In scenario (a) with mi -C m2 <^ m^, this ?7ij-dependence for 6^, (and also for 6) 
is naturally expected, leading to the correlation between two oscillation parameters (e.g., 
between Am-li and sin^ 29^ in the solar neutrino oscillations). For scenario (b) such kind 
of rrij-dependence might not be welcome, since large cancellation between two almost 
degenerate neutrino masses is unavoidable in calculating 6i, and 6. For this reason, we 
like to follow an approach in which the lepton flavor mixing angles do not depend on 
the neutrino masses. This is another strategy for our subsequent introduction about a 
phenomenological model of lepton masses and flavor mixing. 

4.3 Nearly bi-maxiinal mixing 

The dominance of m,- in the hierarchical mass spectrum of charged leptons implies a 
plausible limit in which the mass matrix takes the form 



Ml = c,\ (4.7) 
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with Ci = itIt-. This mass matrix is equivalent to the following matrix with 5'(3)l x S'(3)r 
symmetry or flavor democracy: 



{U 



through the orthogonal transformation I/qMiUq = Mqi, where 

f^o = I 76 75 -^ • (4.9) 

^73 73 73 / 

For either Mi or Mq/, further symmetry breaking terms can be introduced to generate 
masses for muon and electron. Recently the possible signiflcance of the approximate 
democratic mass matrices has been remarked towards understanding fermion masses and 
flavor mixing 



The similar picture is however invalid for the neutrino sector, if three neutrino masses 
are almost degenerate. Provided the neutrino masses persist in an exact degeneracy 



symmetry, the corresponding mass matrix should take the diagonal form |21,22|: 




M, = c J r/2 , (4.10) 



where Ci^ = mo = \mi\ (for i = 1,2,3) measures the mass scale of three neutrinos, and 
1]^ = ±1 (denoting the CP-parity if neutrino masses are of the Majorana type). In the 
case of ?7i = ?72 = ^73, M^y becomes Mqi, which has exact S'(3) sjTumetry and rrii = rriQ: 



Mo. = c J 1 . (4.11) 




By breaking the mass degeneracy of M^, or Mq^ slightly, one may get the realistic neutrino 
masses mi, m2 and ma (at least two of them are different from mo and different from 
each other). 

A purely phenomenological assumption is that the fundamental theory of lepton in- 
teractions might simultaneously accommodate the charged lepton mass matrix Mqi and 
the neutrino mass matrix Mqi, (or more generally. My) in a speciflc basis of flavor space. 
Of course there is no flavor mixing in this symmetry limit. The realistic lepton mass 
spectra and the flavor mixing matrix depend on the explicit introduction of perturbative 
corrections to Mqi and Mq^^. To avoid the dependence of flavor mixing angles on ttij, 
however, only the diagonal perturbations or the special off-diagonal perturbations to Mq^^ 
are allowed. 

We flrst describe an instructive ansatz with the diagonal perturbations to Mqi, (or 
My). The flrst step for symmetry breaking is to introduce small corrections to the (3,3) 
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elements of Mm and Mn,,- The resultant mass matrices read [E3 



i-Ou 



M,i 



Ml, 



1 


1 


1 


1 


1 


1 


1 


1 


1+5, 


1 











1 











l+e^ 



(4.12) 



where \ei\ <^ 1 and \e^\ <^ 1. Now the charged lepton mass matrix ceases to be of rank 
one, and the muon becomes massive (m^ = 2\ei \mr/9 to the leading order of e,). The 
neutrino mass m^ is no more degenerate with mi and 7tt,2 (i.e., |m3 — mo| = ■mol^^l). 
It is easy to see, after the diagonalization of Mn and Mi^, that the second and third 
lepton families have a definite flavor mixing angle 9. We obtain tan^ = — v^ if the small 
correction oiO{m^/mT-) is neglected. Then neutrino oscillations at the atmospheric scale 
arise in v^ <->■ Vj. transitions with the mass-squared difference Amgg = Arrag-^ ^ 2r7io|£;^| 
and the mixing factor sin^ 20 ^ 8/9. Such a result is in good agreement with current 
data. 

The next step is to introduce small perturbations to the (2,2) and (or) (1,1) elements 
of Mil and Mii,, in order to generate the electron mass and to lift the degeneracy between 
mi and m,2. It has been argued in Refs. |lT2|,p!^, in analogy to the quark case, that at 
this step a simple and instructive perturbation to Mu should be of the form that its (1,1) 
and (2,2) elements simultaneously receive small corrections of the same magnitude and of 
the opposite sign. The analogous correction can be introduced to Mi^. Then the lepton 
mass matrices become [113 



M- 



21 



M. 



2v 




(4.13) 



where \bi\ <^ 1 and \by\ <^ 1. One finds me = |5/pm^/(27m^) to the leading order as 
well as mil = ""^0(1 ~ ^u) and m,2 = m,Q{l + 6^). The diagonalization of M21 and M21, 
leads to a full 3x3 flavor mixing matrix, given as f/o in Eq. (9) if small corrections 
of 0{^Jmte/m^j) and 0{m,f^/m,^) are neglected. Then the solar neutrino deficit can be 
interpreted by z/g ^^ '^^i oscillations with the mass-squared difference Arn^i ~ 4?7io|(5i/| 



and the maximal oscillation amplitude [12]. 

If the corrections from non-vanishing muon and electron masses are taken into ac- 
count, the lepton flavor mixing matrix will in general read as V = OiUo, where Oi is an 
orthogonal matrix. Three rotation angles of Oi are functions of the mass ratios mte/m,^ 
and m,n/m,r- Due to the strong hierarchy of the charged lepton mass spectrum |l^ 



I.e., 



a 



(3^ 



m,p 



-^ ^0.0695 



m„ 



m,. 



m^ 



0.0594 , 



(4.14) 
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Oi is expected not to deviate much from the unity matrix. In our specific symmetry- 
breaking case discussed above, we obtain 



h' 




Oi= \ -« 1 - t«^ - iP' -^/5 I (4.15) 

to the next-to-leading order. Note that there is another solution for Oi and it can directly 
be obtained from Eq. (15) with the replacements a -^ —a and f3 -^ —f3. The lepton 
fiavor mixing matrix turns out to be 

V(±) = Uo±{aA-PB)-{a^C-aPD + P^E) , (4.16) 



in which A,---,E are constant matrices and their explicit forms can be found in Ref. |jT3| 
(or read off directly from the product of Oi and Uq). The effects of 0{a'^), 0{aP) and 
Oi^P"^) on neutrino oscillations will be discussed subsequently. 

The flavor mixing matrix obtained above can be parameterized as that in Eq. (2). 
Under CP invariance we are left with only three Euler angles. We then obtain tan^; = 0, 
tan6'jy = 1 and tan 6' = —\/2 in the limit where terms of 0{a) and 0(/9) are neglected. 
Taking small corrections of 0{a) and 0{f3) into account, we arrive at tan^/ = ±a, 
tan6'j, = 1 and tan6' = —\/2 (1 ± 3/5/2), where the "±" signs correspond to V(±) in Eq. 
(16). The full next-to-leading-order results for three mixing angles are found to be 

tan 6^/ = ± a ( 1 =F -/9 

tan^i, = 1 — 3v3 a/5 , 

tan^ = -V2 (l±-[n . (4.17) 

One can see that the rotation angle 9,, only receives a tiny correction from the charged 
lepton masses. 

Following Eq. (1) we take Aml^^ = /S.m\^ and Am^^j^ = Arrigg ~ ^f^\i to accommo- 
date current data on solar and atmospheric neutrino oscillations. Calculating the survival 
probability P{ve —>■ ^e) and the transition probability ^(t'^ -^ Vr) to the next-to-leading 
order, we arrive at [|T3| 

g 

sin^ 26'sun = 1 - -a^ , 

siY?29^,^ = \{lTl3) ■ (4.18) 



This is just the "nearly bi-maximal" mixing pattern, flrst proposed in Ref. ||12[. The 
numerical results for mixing angles in Eqs. (17) and (18) are listed in Table 1. One can 
see that the flavor mixing patterns "V'(+)" and "V(_)" are both consistent with the present 
data on atmospheric neutrino oscillations. For solar neutrino oscillations, our result favor 
the "Just-So" solution. 

The near degeneracy of three neutrino masses assumed above leads to 



7712 — fni 



rris — 1712 



Ayr? 

-—^ ~ 10-^ ("Just -So" solution) . (4.19) 
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Table 4.1: Numerical results for mixing angles 6i, 9,^, 9 and sin^ 26'suiu sin^ 26'atm in the 
"nearly bi- maximal" mixing scenario |T^,p!3|. 



Case 9i 9^ 9 sin^ 26'sun sin^ 26'atm 



"1/f+)" +3.6° 44.4° -57.0° 0.99 0.84 






"1// ^" -4.3° 44.4° -52.2° 0.99 0.94 



This kind of neutrino mass spectrum can account for the hot dark matter of the universe, 
ii TTii K. 2 eV (for i = 1,2,3). The relatively large gap between Am2i and Arrigj (or 
Am|^) has some implications on the forthcoming long-baseline experiments. 

Now we consider the effect of non- vanishing 9i on the survival probability of electron 
neutrinos in a long-baseline (LB) experiment, in which the oscillation associated with the 
mass-squared difference Am^i can be safely neglected due to Am^2 ^ ^^32 ~ ^'^si • It 
is easy to find 

Pii^e - ^e)LB = 1 " ^«^ (1 + 2/5) siu^ f 1.27 ^^) , (4.20) 



The oscillation amplitude amounts to 1.1% (the "V(+)" case) or 1.4% (the "V(_)" case) 
and might be detectable. One can see that the small mixing obtained here lies well 
within the allowed region of sin^ 2^ch from the CHOOZ experiment 0]. The transition 
probability of u^ to z/g in such a long-baseline neutrino experiment reads 



Piu, -. ..),B = y «^ (1 T P) sin^ (l.27 ^^\ . (4.21 



Here the mixing factor is about 0.8% (the "V(+)" case) or 0.9% (the "V(_)" case). The 
proposed K2K experiment is expected to have a sensitivity of sin^ 29 > 10% for u^ ^^ ^^l 
oscillations, while the MINOS experiment could probe values of the mixing as low as 
sin^26' = 1% ||2^. Thus a test of or a constraint on the prediction obtained in Eqs. (20) 



and (21) would be available in such experiments. The K2K experiment can definitely 
measure the survival probability of v^ neutrinos, which reads as 

P{^, - ^.)lb = 1 - ^ (1 + /3) sin^ (l.27 ^^) . (4.22) 

The mixing factor, corresponding to two different perturbative corrections of the magni- 
tude (3 ~ 6%, takes the value 0.84 or 0.94 (see Table 1). It is also worth mentioning that 
the transition probability of v^. -^ z/^, which satisfies the sum rule 

P{l^e -^ Z^e)LB + Pi^e ^ V^.\B + P{l^e ^ l^r)l.B = 1 , (4.23) 
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is smaller (with the mixing factor 8a^/9 ~ 0.4%) and more difficult to detect. 

In the discussions made above the type of neutrinos was not specified. If they are 
of the Majorana type, then their masses have to fulffil the bound (m,^) < 0.45 eV (at 
the 90% confidence level) from the neutrinoless /5/3-decay p5|, where {rrti,) is an effective 



mass factor. The magnitude of {mj) depends both on the neutrino masses mi and the 
flavor mixing elements Vef {m^) = ^ iVi^^i^^) fo^ ^ = 1,2 and 3, where rj^ is the CP 
parity of the i-th Majorana neutrino fleld as illustrated in Eq. (10). Taking rj-^ = —%, 
one flnds that the flrst and second terms in {rriu) cancel essentially with each other due 
to the high degeneracy between mi and m2 as well as the approximate equality between 
V^i and Vg|. The magnitude of {m^} turns out to be 

(m^) ^ —= a rrii . (4.24) 

v3 

We then arrive at rrii ^ 12.4 {m^) < 5.6 eV, a bound well within our original expectation 
rrii r^ 2 eV, which can account for the hot dark matter of the universe. 

An interesting theoretical understanding of our "nearly bi-maximal" mixing ansatz 
has recently been made by Mohapatra and Nussinov on the basis of a left-right symmetric 



extension of the standard model with 5(3) and Z{A) x Z(3) x Z(2) symmetries ||26[ . 



4.4 The "small versus large" mixing scenario 

Now we turn to a different symmetry-breaking scenario for the charged lepton mass 
matrix Mqi and the neutrino mass matrix Mq^. Keeping the diagonal symmetry-breaking 
chain Mq/ — > Mn -^ M21 unchanged, we introduce the off-diagonal perturbations to 
Mo^y. To ensure the "maximal calculability" for the neutrino mass matrix, we require 
a special form of the perturbative matrix: it only has two unknown parameters (to 
break the mass degeneracy of Mov) and can be diagonalized by a constant orthogonal 
transformation (independent of the neutrino masses). Then we are left with only three 
perturbative patterns satisfying these strong requirements, as hsted in Table 2. They can 
be diagonalized by three Euler rotation matrices Rij with the rotation angles dij = 45° 



in the (1,2), (2,3) and (3,1) planes respectively (see also Table 2). In Ref. |]T4[ pattern 
(I) was first proposed and discussed. For each pattern the resultant flavor mixing matrix 
reads as V = OiUoRij, where Uq and Oi have been given in Eqs. (9) and (15). To 
the leading order [a = (3 = 0), V remains a constant matrix as shown in Table 2. We 
calculate the transition probability for i/^ -^ Vr and flnd the mixing factors to be 8/9, 
2/9 and 2/9, corresponding to patterns (I), (II) and (HI). Therefore only pattern (I) can 
survive when confronting the atmospheric neutrino data. 

Let us discuss the consequences of pattern (I) in some detail. Note that rj^ = ri2 = % = 
1 has deflnitely been taken |0. Three neutrino masses are given by mi = mQ{l — e^), 
^2 = ''^0(1 + £u) and ma = mo(l -|- 6u), respectively. To the next-to-leading order the 
flavor mixing matrix V has two slightly different forms, because of two possible solutions 
of Oi- Note that in Ref. 0] only one solution was numerically given. Here we obtain the 
analytical results of three mixing angles for both possibilities: 

tan 61/ = ±a ( 1 =F -/3 
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Table 4.2: Three perturbative patterns for Mqi^ and their consequences on V in the leading 
order approximation (a = = 0) . 



Pattern Perturbation 



Rotation i?, 



Flavor mixing V 



(I) 








1 





(11) 





1 





1 



V2 y/2 

1 1 





(111) 




1 

2 

1 

2v^ 
2 

V6 



1 

V2 
1 

1 
V3 




tan^^ = —a(5 [^Ti^P 



(4.25) 



Comparing these results with those obtained in Eq. (17), one can see that only the value 
of tan^^ changes. Accordingly the mixing factors of solar and atmospheric neutrino 
oscillations read: 



sm 



29. 



sm 



'29. 



9 



;iT/3) 



(4.26) 



to the next-to-leading order. This is just the so-called "small versus large" mixing 
scenario. The numerical results for mixing angles in Eqs. (25) and (26) are listed in 
Table 3. We observe that this ansatz favors the small-angle MSW solution to the solar 
neutrino problem. Its consequence on the atmospheric neutrino oscillations is the same 
as that obtained in the last section. 

It is remarkable that this ansatz also has the same implications on the long-baseline 
neutrino experiments at the atmospheric scale. This can be seen clearly from Eqs. (5) and 
(6), in which the u^^ — > z/g and u^. —>■ Vj. transition probabilities are absolutely independent 
of the mixing angle 9^. Therefore the results obtained already in Eqs. (20) - (22) remain 
valid in the "small versus large" mixing scenario (for a numerical study of these long- 
baseline transitions, see Ref. 12 
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Table 4.3: Numerical results for mixing angles 6i, 9^, 9 and sin^ 20suiu sin^ 20atm to the 
next-to-leading order in the "small versus large" mixing scenario. 



Case 



9, 



9u 



9 



sin^ 26's„r, sin^ 29» 



"TA . " 



'V( 



(+) 



"1/r- ■ " 



'(-) 



+3.6° -0.60° -57.0° 8.0 x lO'^ 
-4.3° -0.63° -52.2° 4.9 x lO^^ 



0.84 
0.94 



If neutrinos are of the Majorana type, then the smallness of both \9i\ and \9j^\ in this 
ansatz implies that the effective mass factor of the neutrinoless /?/3-decay (i.e., {frty)) is 
dominated by mi. A strong constraint turns out to be mi ~ {m,^) < 0.45 eV. In addition, 
the near degeneracy of three neutrino masses assumed above leads to 



7712 — "^1 



f^3 ~ 1^2 



Amt 



^'^atm 



10" 



(4.27) 



Therefore the sum of three neutrino masses has an upper bound of 1.4 eV, which seems 
difficult (if not impossible) to account for the hot dark matter of the universe. 

4.5 Bi-maximal mixing 

Finally let us give some comments on the "exactly bi-maximal" mixing scenario of three 



neutrinos, which is recently proposed by Barger et al |]T5[. The relevant flavor mixing 
matrix, similar to Uq in Eq. (9), reads as follows: 



v 



i 


i 





V2 


V2 


1 


1 


1 


2 


2 


v^ 


1 


1 


1 


2 


2 


^ 



(4.28) 



This flavor mixing pattern is independent of any lepton mass and leads exactly to 
sin^ 26'atm = sin^ 26'sun = 1 for neutrino oscillations. Therefore it favors the "Just-So" 
solution of the solar neutrino problem. We flnd that V can be derived from the following 
charged lepton and neutrino mass matrices [O]: 



m; 



M'. 















(4.29) 



where \5[J\ <^ 1 and \s'iJ\ <^ 1. In comparison with the democratic mass matrix Mq/ 
given in Eq. (8), which is invariant under the S'(3)l x S'(3)r transformation, the matrix 
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M/ in the limit S'l = e'l = only has the >S'(2)l x 5(2)^ symmetry. However Ml in the limit 
SI = el = takes the same form as Mqi, in Eq. (11), which displays the S{3) symmetry. 
The off-diagonal perturbation of M/ allows the masses of three charged leptons to be 
hierarchical: 



c' 



{m, ,m,,m^} = ^ m , \e[\ , 2 + e[} . (4.30) 

We get c'l = m^ + rriT- ^ 1.88 GeV, \e'i\ = 2m^/(m^ + m^) ^ 0.11 and \6'i\ = 2me/(m^ + 
rrir) ~ 5.4 x 10~^. The off-diagonal perturbation of M'^ makes three neutrino masses 
non-degenerate: 

{mi , 1712 , "^3} = 4 {1 + 4 5 1 - 4 5 1 + 4} • (4-31) 

Taking Am'^^^ = Am2i and Aml^^ = Amgg ~ Arrigj^ for solar and atmospheric neutrino 
oscillations, respectively, we then arrive at the same result as that obtained in Eq. (19). 
The diagonalization of M/ and M^ leads straightforwardly to the flavor mixing matrix V. 
In Ref. [|l^ a different neutrino mass matrix has reversely been derived from the given 
V in a flavor basis that the charged lepton mass matrix is diagonal. The emergence of 
the "exactly bi-maximal" mixing pattern from M/ and M'^ in Eq. (29) is, in our point 
of view, similar to that of the "nearly bi-maximal" mixing pattern from M21 and M2U in 
Eq. (13). 

Note that three mixing angles of V are given as ^^ = 0, ^j, = 45° and 6 = —45°. 
The vanishing 61 leads to vanishing probabilities for u^ — > i^^ and u^ — > Uj. transitions in 
the long-baseline experiments at the atmospheric scale. This feature makes the "exactly 
bi-maximal" mixing ansatz experimentally distinguishable from the "nearly bi-maximal" 
mixing ansatz. 

If three neutrinos are of the Majorana type, then the near degeneracy of their masses 
implies that (m^) ~ rrii for the neutrinoless /3/9-decay (note that ?7i = ''72 = V3 has been 
taken in M^). Therefore nii < 0.45 eV and the sum of three neutrino masses seems 
insufficient to account for the hot dark matter of the universe. 

4.6 Conclusions 

In summary, we have discussed a simple phenomenological model of lepton mass gener- 
ation and flavor mixing on the basis of flavor democracy for charged leptons and mass 
degeneracy for neutrinos. Three instructive symmetry-breaking scenarios of this model 
lead to three interesting mixing patterns of three neutrinos, i.e., the "nearly bi-maximal" 
mixing, "small versus large" mixing and "exactly bi-maximal" mixing. 

The "nearly bi-maximal" mixing and "exactly bi-maximal" mixing ansatze favor the 
"Just- So" solution of the solar neutrino problem |l[]. It is possible to distinguish between 
these two scenarios in the future long-baseline neutrino experiments. On the other hand, 
the "small versus large" mixing ansatz favors the small-angle MSW solution to the solar 
neutrino problem. Whether the solar neutrino deflcit is attributed to the "Just-So" oscil- 
lation or the matter-enhanced oscillation can flnally be clarifled by the Super-Kamiokande 
experiment and other neutrino experiments (e.g., the SNO experiment |]2S|) under way. 

As we have emphasized before, our purely phenomenological approach can only be 
"justifled" by its consequences. If it is experimentally favored, one may get useful hints 
towards deeper (theoretical) understanding of the phenomena of neutrino masses and 
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flavor mixing. Recently a lot of attention has been paid to the textures of lepton mass 
matrices and their implications on neutrino physics ||29| , in which some similar symmetry 



arguments were made. Further attempts in this direction are no doubt desirable. 

4.7 Discussion 

Giunti: Which is the difference between your approach and that of Kang et al.l 

Xing: The approaches of Kang et al. (see Ref. |29[) and ours start essentially from a 
similar point of view. The main difference lies in how to introduce perturbative 
corrections. In their case the resulting flavor mixing depends on neutrino masses, 
while in our case the independence of flavor mixing on neutrino masses is required 
from the beginning. The results of both approaches are consistent with current 
data. 
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5 Field theoretic treatment of source lifetime and 
wave-packet effects in neutrino oscillations — Sub- 
hendra Mohanty 

In the wave-packet treatments of neutrino oscillations the initial neutrino wave-packet is 
assumed to be a gaussian with spread ax in the initial position, and at the uncertainty in 
the initial time of production of the particle. For accelerator neutrinos which are produced 
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from the decay of pions or muons, it is assumed that at = v^ x Tsource- We show below 
that the gaussian wave-packet treatment for the case of long lifetime neutrino sources 
— like muons in the LSND experiment — does not reproduce the quantum mechanics 
formula for neutrino oscillations. 

The propagation amplitude of a Gaussian wave-packet turns out to be [2], 



ir(X, T; ma) = exp{-iEaT + iP^ ■ (X)} exp 



(X-v.T)^ 



(5.i: 



4(^^ + v^af) j ' 

where X = Xf — Xi, T = t/ — tj and (a^ + f^cr^)^/^ = o" is the width of the wave-packet. 
The probability of flavor oscillation as a function of space-time is 

2 



P{a^ P;X,T) 



3 
a=l 



U^aK{ma;X,T)U* 



(5.2) 



The probability as a function of distance is given by the time average of P{a -^ /3; X, T): 



P{a^ P;X) 



dT 



Y,U.aK{ma;X,T)U:^, 



a=\ 



The interference term is 
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where the exponential damping factor A is given by 



A = i^Ef — + 



a 

2^^ 



AP 



X^ 
4^ 



(5.3) 



(5.4) 



(5.5) 



The specific form of the phase difference EAE — P ■ AP which appears in a covariant 
calculation has interesting applications. One can show that for this form of phase dif- 
ference there are no EPR-type associated particle oscillations [2]. Consider a pion decay 
n ^ fiu in the rest frame of the pion. Using the conservation laws 



E„ 



E„ 



ml — m^ -I- ml 

2m,T^ 
ml + ml — ml 



P 



2m7r 

IP I = IP 



[{ml - (m^ + m^,f){ml - (m^ - m^f)]^^"^ 



2m,Tr 



the neutrino phase differences are given by 

Ami 
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2AP. 
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Ami /rn X 

( El, — m^ ) 



X Am^ 
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The phase differences between different muon states are given by 



AE„ 



AP„ 



Ami 

2r7i7r 



A0^ = {E^AE^ - P,, • AP, 



X 



0. 



This shows that there are no EPR-type secondary particle oscillations. 
Going back to the neutrino oscillation term ( p.4|) , we consider the case 



- _ / 2 I 2 2\l/2 

When the source-detector distance X is such that 

2aP 



VT . 



X -C Xcoh 



AP ' 



(Aw? r\ ' 



we have 
The resulting wave-packet formula for neutrino oscillation is 



(5.6) 



(5.7) 



(5.^ 



P[v^ ^ V(,; X) = - sin^ 2^ a - cos ( 



exp 



l.TQAmV 

E 



(5.9) 



In the LSND experiment (neutrinos from muon decay) the source lifetime is r^ = 
658.6 m and the exponential factor is significant. 

In other accelerator experiments with neutrinos from pion or kaon decay the source 
lifetimes is ~ 7 m and the exponential factor is negligible. 

From Figure p.l| it is clear that the standard gaussian wave-packet treatments does 
not reproduce the quantum mechanical formula when the source of the neutrinos has a 
large lifetime. 

Why is the neutrino wave-packet treatment inconsistent? 

A spatial width a^ ~ 659 m of the neutrinos implies that the spread of neutrino 
momentum is ap = l/2cr^ ~ 10~^^ MeV. This is experimentally unreasonable as the 
spread of the neutrino momentum must be at least as large as the momentum spread of 
the source (pion, muon) wave packet, which is expected to be ~ 1 — 10 MeV. 

To incorporate the source lifetime we must include the Breit-Wigner propagator of 
the source in the Feynman diagram [3]. 

We consider the LSND experiment with the process 



/^ 



+ i^e + i^fj. ^^' Ve + P ^ n + e^ 



(5.10) 



as a model for our investigation of the source- lifetime dependence of neutrino oscillations. 
The amplitude of the process {11^)3 + {p)d -^ (e"'"i^e)s + (n -|- e"*")/) is expressed as the 
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The cross section obtained by taking the modulus square of ( 5.11 
oscillation probability [4]: 
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• In the limit of a small source- lifetime (F ^ Ai?), which is valid for example for 
neutrinos form Z decay, the lifetime dependence drops out and we get back the 
standard oscillation formula. 

• The most interesting application of (|5.12|) is in the regime F < !\E. In this regime, 
which includes LSND, the Am^ probed by ( |5.12| ) is lower than what one obtains 
from the standard oscillation formula. 

• In earlier treatments os and ar, were taken to be the momentum spread of the 
neutrino wave-packet, that is a quantity over which there is no experimental control. 
In this analysis as and gd are the spread of the muon and proton momenta which 
are better known experimentally. 

5.1 Discussion 

Giunti: I think that the sensitivity to IS.w? of a neutrino oscillation experiment is de- 
termined only by the argument Aim?X/E of the oscillatory term. If Pmax is the 
maximum value of the oscillation probability measurable in a given experiment, 
the minimum value of Am^ that can be probed is Atti^j^ ~ V-Pmax-E/X (with 
appropriate averages over energy and distance). 

A long lifetime of the source means that the coherence length is very large, but it 
cannot have any effect on Am^jj^. Indeed, a long lifetime of the source implies that 
the spatial width of the neutrino wave packet is very large and the standard plane- 
wave treatment is applicable. A correct wave-packet treatment must reproduce this 
result. 

In the case of the LSND experiment the anti-muon does not decay exactly at rest, 
because it is in a medium and its kinetic energy must be of the order of the thermal 
energy of the medium. The corresponding velocity is f ^ ~ 3 x 10~^. The mean free 
path of the anti-muon in the medium is of the order of the inter-nuclear distance: 
£^ ~ 10~^ cm. Therefore, the coherent emission of the neutrino wave is interrupted 
every (!-f_i/v^ ~ 10^^^ s, which is much smaller than the muon lifetime, r^ ~ 2 x 10^^ s. 

Furthermore, the positron emitted in the anti-muon decay /i"*" -^ e~^ + z/g + z/^ is 
relativistic {v^ — 1) and is annihilated as soon as it interacts with an electron in 
the medium. Approximating its mean free path with the inter-nuclear distance, 
4 ~ 10~^ cm, one can see that the coherence of the //"""-decay process is interrupted 
after ie/ve ~ 3 x 10~^^s. This is the dominant effect in the determination of 
the spatial width of the neutrino wave packet, which results to be of the order of 
10~^ cm, much smaller than the ~ 600 m obtained from the muon lifetime. 



'■1 



Mohanty: I agree that in LSND the muon decay does not happen in vacuum and the 
correct value for at should be the mean free path of the muon. 

My analysis is not specific to LSND but to the general case of neutrinos produced 
from isolated long lived particles. Then, from the oscillation formula (5^) one can 
see that when there is an incoherent superposition (for example when the detector 
distance is much larger than the coherence length) we have P = \ sin^ 29. The plot 
of this probability in the /S.rn? vs sin^ 26 plane gives a vertical line at | sin^ 20 for 

52 



all values of Am^ for which the decoherence condition is satisfied. So the general 
result is that the incoherent superposition formula rules out a larger portion of the 
Am^ vs sin^ 29 graph than the coherent superposition formula. Therefore, it is of 
practical interest to examine various ways in which we can obtain an incoherent 
superposition. 

The second point I wanted to make is that the assumption of a gaussian wave-packet 
with width equal to the lifetime, which is the standard assumption in literature (see 
Ref. and references therein), does not reproduce the quantum mechanics formula. 

Lipkin: There is one problem with the "field theoretic treatment" of the source lifetime. 
The problem is that the neutrino is not emitted alone from the source; there is also 
another change in the environment. If we are considering a long-lived beta decay of 
a nucleus bound in an atom, the nuclear lifetime is irrelevant for neutrino coherence 
because the nucleus is interacting with the atom, and the atom knows when the 
charge of the nucleus has changed and an electron or positron has been emitted 
together with the neutrino. 

The point has been repeatedly made by Leo Stodolsky that the proper formalism 
to treat neutrino oscillations is not field theory but the density matrix, because 
only in this way the unavoidable interactions with the environment can be taken 
into account. Leo also points out that the length in time of the wave packet is 
irrelevant 0. 
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6 The Nucleosynthesis Limit on N^ — Subir Sarkar 

Hoyle and Tayler |jl| as well as Peebles had emphasized many years ago that new 
types of neutrinos (beyond the z/g and u^ then known) would boost the relativistic energy 
density hence the expansion rateQ during big bang nucleosynthesis (BBN), thus increas- 
ing the yield of ^He. Shvartsman [Q noted that new superweakly interacting particles 
would have a similar effect. Subsequently this argument was refined quantitatively by 
Steigman, Schramm and collaborators [ffl . In the pre-LEP era when the laboratory bound 



^In the radiation-dominated era, H = -^/SttGn/o/S, with p — ^g^.T'^ where g^, counts the relativistic 
degrees of freedom. 
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on the number of neutrino species was not very restrictive 0, the BBN constraint al- 
ready indicated that at most one new family was allowed 0, albeit with rather uncertain 
systematics ^. Although LEP now finds N^, = 2.991 ±0.016 p|, the cosmological bound 
is still important since it is sensitive to any new light particle, not just SU{2)l doublet 
neutrinos, so is a particularly valuable probe of new physics 0.0 

The primordial mass fraction yp(^He) increases as ~ 0.012AA^j^ but it also increases 
logarithmically with the nucleon density (usually parameterized as r^ = n-^/n^ = 2.728 x 
lO^^Q-^h"^) . Thus to obtain a bound on N^, requires an upper limit on Yp and a lower 
limit on t]. The latter is poorly determined from direct observations of luminous matter 
so must be derived from the abundances of the other synthesized light elements, D, ^He 
and ^Li, which are power-law functions of t]. The complication is that these abundances 
are substantially altered in a non-trivial manner during the chemical evolution of the 
galaxy, unlike yp(^He) which just increases by a few percent due to stellar production. 
(This can be tagged via the correlated abundance of oxygen and nitrogen which are made 
only in stars.) 

Even so, some cosmologists have used chemical evolution arguments to limit the pri- 
mordial abundances of D and ^He and thus derived increasingly severe bounds on N^, ||10[ , 



culminating in a recent one below 3 [rT|! However a more conservative view [0 is that 
there is no crisis with BBN if we recognize that such arguments are rather dubious and 



consider only (izrect measurements |I3| of light element abundances, as shown in Figure 1. 
The '^He mass fraction is obtained from observations of metal-poor blue compact galax- 



ies by linear extrapolation to zero nitrogen/oxygen abundance [0; the upper limit is 
reliable, the lower one less so. In particular, older measurements had suggested a value 
smaller by about 4% [llH]. At present there are two conflicting measurements of the D 



abundance in quasar absorption systems |T6|,IT^; the higher value |T6| is interpreted as an 



upper limit. Also shown is the abundance in the interstellar medium ||18[ which provides 
a reliable lower limit. The ^Li abundance as measured in the hottest, most metal-poor 



halo stars [IS] as well as in disk stars |20] is shown and interpreted as providing, respec- 
tively, reliable lower and upper limits on its primordial value. Given these uncertainties, 
standard BBN is consistent with observations for 77 ?^ 2 — 9 x 10^^°. Adopting the reli- 
able limits, yp(^He) < 0.25, D /H > 1.1 x 10"^ and '^Li/H < 1.5 x 10"^ and taking into 
account uncertainties in nuclear cross-sections and the neutron lifetime by Monte Carlo, 



we obtain 12 



iV^^- = 3.75 + 78 (Fp™^" - 0.240), (6.1) 

i.e. up to 1.5 additional (equivalent) neutrino species are allowed for rj at its lowest 
allowed value. It is clear that the restrictions on new physics are less severe than had 
been reported previously |1T0|. 



Other workers have applied Bayesian likelihood methods to obtain Ni, < 4 — 5 pT 
In order to enable extraction of the best-fit N^/t] values as the observational situation 
improves (and uncertainties in the input nuclear cross-sections decrease), we have de- 
veloped a simple method for determining the (correlated) uncertainties of the expected 
abundances [^|. Essentially we compute the full covariance error matrix, having checked 



the linearity of error propagation, and then perform a simple x fit to the observational 



■^The energy density of new light fcrniions i is equivalent to an effective number A Niy = 
^j(gi/2)(ri/ri,)^ of additional doublet neutrinos, where Ti/T^, can be calculated from considerations of 
their (earher) decouphng [||. 
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Figure 6.1: Predicted light element abundances for the Standard Model versus the nucleon- 



to-photon ratio |12|. The 95% c.l. limits determined by Monte Carlo reflect the uncertainties 
in input nuclear cross-sections and the neutron lifetime. Rectangles indicate observational 
determinations and associated '95% c.l.' bounds. 



data. This requires us to input actual measurements, not limits as above so we consider 
both the reported (and mutually incompatible) values for the deuterium abundance: 



D/H 



3.40 ± 0.25 X 10"^ (low) 
1.9 ±0.5x10"^ (high), 



as well as for the helium mass fraction: 



yp(^He) = 0.234 ±0.0054 (1 



ow 



0.244 ±0.0054 (high). 



(6.2) 
(6.3) 



(6.4) 
(6.5) 



In Figure 2 we show the likelihood contours in the Ny-rj plane for the four possible 
combinations of the above measurements. (The results are relatively insensitive to the ^Li 
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abundance.) We see that both the high D/low '^He and the low D/high ^He combinations 
are consistent with Ni, = 3 and indeed the x^ ior the fits is excellent |2^ . However the low 
D/high ^He combination suggests that N^ ~ 2.3, while the high D/low ^He combination 
suggests that Ni, ~ 3.8! If either of these possibilities are indeed substantiated by further 
measurements, then this would indicate a departure from the Standard Model, e.g. the 
presence of a sterile neutrino which mixes with the left-handed doublets. Interestingly 
enough, this can both raise and lower the effective value of N^, [ ^3| | . 

A clear discrimination between these possibilities will be provided when forthcoming 
precision measurements of small-scale angular anisotropy in the cosmic microwave back- 
ground ||2^ provide an independent measure of the nucleon density t] to an accuracy of 



a few per cent [P5|] . 
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7 Are the results of LSND and KARMEN 2 com- 
patible? — Carlo Giunti 

The KARMEN collaboration has reported recently [0] a null result of the KARMEN 2 
experiment searching for neutrino oscillations in the channel u^ ^ Ue- This result is very 
interesting because the KARMEN neutrino oscillation experiment [@,|^ is sensitive to 
the same region in the plane of the neutrino mixing parameters sin^ 29 and Am^ as the 
LSND experiment @] whose results provide an evidence in favor of z/^ —>■ Ve and v^ — >■ v^ 
oscillations. Hence, the statistical interpretation of the KARMEN 2 null result (as well as 
that of the LSND result) is crucial in order to obtain an indication on the compatibility 
or incompatibility of the different results of the two experiments. 

So far the KARMEN 2 experiment measured no events |l|, with an expected back- 
ground of 2.88 ± 0.13 events. This null result has been analyzed with the following 
statistical methods: 

Bayesian Approach. This method is accepted by the Particle Data Group [^, H and 
has been used by the KARMEN collaboration [|I| . The resulting upper limit for the 
mean /i of neutrino oscillation events is 2.3 and the corresponding exclusion curve 
is shown in Fig. |7.1| (the solid curve passing through the filled squares). 

Unified Approach. This frequentist method has been proposed recently by Feldman 
and Cousins |^. It is very attractive because it allows to construct classical confi- 
dence belts with the correct coverage that "unify the treatment of upper confidence 
limits for null results and two-sided confidence intervals for non-null results" [Q. 
On the other hand, the Unified Approach has the undesirable feature that when 
the number of observed events is smaller than the expected background, the upper 
limits for the mean /x of true neutrino oscillation events decreases rapidly when the 
background increases. From a physical point of view this is rather disturbing, be- 
cause a stringent upper hound for fj, obtained by an experiment which has observed 
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a number of events significantly smaller than the expected background is not due to 
the fact that the experiment is very sensitive to small values of fi, but to the fact 
that less background events than expected have been observed. 

This is the case of the null result of the KARMEN 2 experiment, from which the 
Unified Approach yields an upper 90% CL confidence limit of 1.1 events for the 
mean /z of neutrino oscillation events. The corresponding exclusion curve in the 
plane of the neutrino mixing parameters sin^ 29 and Ani^ is shown in Fig. ^]l] (the 
solid curve passing through the filled circles) and one can see that it is significantly 
more stringent than the exclusion curve obtained with the Bayesian Approach (the 
solid curve passing through the filled squares). The strictness of the Unified Ap- 
proach exclusion curve is due to the non-observation of the expected background 
events and not to the sensitivity of the experiment (see the discussion in Ref. Q 
and the sensitivity of the KARMEN experiment presented in Ref. 0). This is 
clearly an undesirable result from a physical point of view, because the statistical 
interpretation of the data produces an exaggeratedly stringent result that could lead 
to incorrect physical conclusions. 

The discrepancy between the Bayesian Approach exclusion curve and the Unified 
Approach exclusion curve is worrying for a physicist, because the Bayesian Approach 
exclusion curve is compatible with a large part of the LSND-allowed region (the shadowed 
area in Fig. [7.1|) , whereas the Unified Approach exclusion curve excludes almost all the 
LSND allowed region. 

In view of the uncertainty of the physical meaning of the KARMEN 2 null result in- 
duced by the significant difference between the exclusion curves obtained with the Unified 
Approach on one hand and with the Bayesian Approach on the other hand, it is inter- 
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esting to explore other possibilities for the statistical interpretation of the KARMEN 2 
null result. 

In this report I present the outcome of the statistical analysis of the null result of 
the KARMEN 2 experiment with the New Ordering Approach that has been proposed 
in Ref. 0. This approach is is based on a new ordering principle for the construction 
of a classical frequentist confidence belt that has all the desirable properties of the one 
calculated with the Unified Approach and in addition minimizes the effect on the resulting 
confidence intervals of the observation of less background events than expected. Hence, 
it is appropriate for the statistical interpretation of the null result of the KARMEN 2 
experiment. The resulting upper limit for the mean /x of true neutrino oscillation events 
in the KARMEN 2 experiment is 1.9 and the corresponding exclusion curve [^ is shown 
in Fig. ^]2] (the solid curve passing through the filled triangles). 

The exclusion curve obtained with the New Ordering Approach lies close to the 
bayesian exclusion curve and tends to support the compatibility of the KARMEN 2 
and LSND results. This is a desirable achievement. Furthermore, it is important to 
emphasize that the New Ordering Approach gives a correct frequentist coverage as the 
Unified Approach. 

Hence, the New Ordering Approach has solved the apparent conflict between the 
frequentist and bayesian statistical interpretation of the null result of the KARMEN 2 
experiment: by choosing an appropriate ordering principle in the construction of the 
confidence belt, the exclusion curve obtained with the frequentist method is in reasonable 
agreement with the one obtained with the Bayesian Approach. 

The different impact of the non observation of the expected background events on 
the exclusion curves obtained with the Unified Approach and with the New Ordering 
Approach is illustrated in Figs. U^ and |7.4|. In these figures the exclusion curves obtained 



60 



from the null result of the KARMEN 2 experiment (n = 0, where n is the number of 
observed events) are compared with the exclusion curves that would be obtained by an 
imaginary experiment similar to the KARMEN 2 experiment but with a sensitivity twice 
that of the KARMEN 2 experiment, which observe two events (n = 2), in good agreement 
with the expected mean background h = 2.88 ± 0.13. The results of both experiments 
(real and imaginary) are compatible with the absence of neutrino oscillations and imply 
upper bounds for the neutrino oscillation parameters. 

One can see from Fig. [7.3| that the imaginary experiment with twice sensitivity would 
yield a Unified Approach exclusion curve (the solid curve passing through the empty 
circles) which is significantly worse than that obtained from the null result of the KAR- 
MEN 2 experiment (the solid curve passing through the filled circles) . On the other hand 
Fig. \7.4\ shows that the imaginary experiment with twice sensitivity would yield a New 
Ordering exclusion curve (the solid curve passing through the empty triangles) which is 
practically equivalent to that obtained from the null result of the KARMEN 2 experiment 
(the solid curve passing through the filled triangles). It is clear that the impact of the non 
observation of the expected background events on the resulting exclusion curve cannot 
be eliminated completely, but can be minimized with an appropriate construction of the 
frequentist confidence belt, as that resulting from the New Ordering Principle. 

I conclude this report with the following remarks: 

1. The statistical interpretation of the null result of the KARMEN 2 neutrino oscil- 
lation experiment is rather problematic because no events were observed with a 
mean expected background of 2.88 ±0.13 events |l|]. The exclusion curves obtained 
with the Bayesian Approach and with the Unified Approach [0 are significantly dif- 
ferent and yield contradicting indications on the compatibility of the KARMEN 2 
result with the neutrino oscillation signal measured in the LSND experiment |^ 
(see Fig. O). 



2. The analysis of the KARMEN 2 null result with the New Ordering Approach [^, 
which is a frequentist method with correct coverage as the Unified Approach, yields 
an exclusion curve close to the one obtained with the Bayesian Approach (see 
Fig. |7.2| ). In this way, the undesirable discrepancy between frequentist and bayesian 
interpretations of the KARMEN 2 null result is removed. 

3. Taking into account the error of the expected mean background b = 2.88 ± 0.13 
in the KARMEN 2 experiment, even if it is wrong by an order of magnitude, does 
not help in solving the problem of the statistical interpretation of the result of this 
experiment, because the resulting exclusion curves are practically equivalent to the 
ones obtained assuming no error for the expected mean background [ |T0[] . 

4. An extreme attitude^ is to ignore the calculated mean expected background and to 
assume that the background is unknown. This approach gives ultra-conservative 
exclusion curves [|1^], which in the case of the KARMEN 2 experiment tend to 



*Let me emphasize that this attitude is purely speculative and does not have any justification from the 
experimental point of view. The background in the KARMEN 2 experiment (except for the component 
due to the intrinsic 9e contamination of the beam, that is expected to contribute with 0.56 ±0.09 events) 
is measured on-line with high precision in parallel with the neutrino oscillation search jl^, n2| . 
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support the exclusion curve obtained with the Bayesian Approach. Obtaining ultra- 
conservative exclusion curves is generally not desirable, but could be considered to 
be a safe choice in controversial cases as that of the KARMEN 2 experiment. 

5. The sensitivity of an experiment has been defined by Feldman and Cousins as 
"the average upper limit that would be obtained by an ensemble of experiments 
with the expected background and no true signal" 0. The 90% CL sensitivity 
curves obtained from the null result of the KARMEN 2 experiment with the Uni- 
fied Approach and with the New Ordering Approach are shown in Figs. ^!T] and 
[77^ (the sohd curves passing through the empty circles and triangles, respectively). 
In Ref. [0 it has been shown that a sensitivity curve cannot be considered as an 
exclusion curve. Since the sensitivity curve of a neutrino oscillation experiment can 
be calculated before doing the experiment, its usefulness lies in the possibility to 
plan future experiments in order to cover approximately the region of interest in 
the plane of the neutrino mixing parameters sin^ 29 and Am^. However, a large dis- 
crepancy between the exclusion and sensitivity curves is a signal that the exclusion 
curve results form a rather improbable experimental result and one must be very 
cautious in formulating physical conclusions on the basis of the exclusion curve. It 
can be seen from Figs. |7.1| and |7.2| that from this point of view the new ordering 
approach is safer than the unified approach. 

6. The direct comparison of exclusion curves and allowed regions obtained with dif- 
ferent statistical methods does not have a precise statistical significance. Hence, 
such a comparison cannot be used to combine the results of different experiments 
or to infer with some known confidence level a contradiction between the results of 
different experiments when the comparison is done on the border of the exclusion 
curves and of the allowed regions. Hence, the comparison of the KARMEN 2 ex- 
clusion curves and the LSND-allowed region, which was obtained with a different 
statistical analysis (see Ref. |Q), must be done with great caution. 

7. It is clear that the null result of KARMEN 2 does not favor the credibility of the 
positive result of the LSND experiment, but I think that it is still premature to 
claim a contradiction between the two experiments. 

8. In the future, if the KARMEN 2 experiment will continue to observe no neutrino 
oscillations, it will be possible to claim a contradiction between the results of the 
two experiments only when the exclusion curves obtained from the result of the 
KARMEN 2 experiment with different statistical methods will produce similar re- 
sults and will lie well on the left of the region allowed by the results of the LSND 
experiment. 

7.1 Discussion 

Sarkar: I think that the Unified Approach exclusion curve is fine, as long as it is inter- 
preted in the correct way. 

Giunti: Yes, I agree. However, the majority of physicists are not expert in the subtleties 
of statistics and tend to believe in the exclusion of the parameter region on the right 
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of an exclusion curve. This is statistically incorrect and very misleading in some 
cases, as that of the KARMEN 2 experiment. Therefore, I think that it is desirable 
to present the experimental results in a way that leaves less space to improper 
interpretations. I think that this can be done using the New Ordering Principle, 
which is equivalent to the Unified Approach from the statistical point of view: both 
of them give an allowed range of the mixing parameters that belongs to a set of 
allowed ranges that could be obtained with an ensemble of experiments identical 
to the KARMEN 2 experiment and cover the true values of the mixing parameters 
with probability 0.90 (for 90% CL exclusion curves). 
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